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HARTREK-rOCK  THEORY  IH  NUCLEAR  PHYSICS 


D>  Go^y*  and  P.  L<  Llona** 

I.  Introduction 

This  papar  is  devoted  to  a  general  presentation  of  Hartree-Pock  equations  and  related 
questions,  and  we  will  be  mainly  Interested  In  the  application  of  Rartree-Fock  method  to 
Nuclear  Physics. 

As  It  Is  well-known,  the  Rartree-Pock  method  was  Introduced  by  D.  Hartree  [23],  V. 

Pock  [17]  and  J.  C.  Slater  [45]  to  approximate  the  ground  state  (and  Its  energy)  of  general 
N-body  problems  In  Quantum  Physics.  And  the  main  application  of  this  method  was.  In  Atomic 
Physics,  the  study  of  Coulomb  systems  (atoms  and  molecules)  with  the  purely  Coulomb 
Hamiltonian  of  electrons  Interacting  with  static  nuclell. 

In  Nuclear  Physics,  the  use  of  Hartree-Pock  methods  to  compute  the  ground  state  of 
nuclell  is  quite  recent  (see  for  example  the  review  papers  by  H.  Bethe  [7],  J.  w.  Negele 
[41],  [42],  P.  Quentin  and  H.  Plocard  [44]  and  the  references  therein);  among  other 
reasons,  this  delay  was  due  to  the  lack  of  understanding  of  strong  Interaction  and  thus  to 
the  difficulty  of  deriving  realistic  Hamiltonians  to  describe  the  Interaction  of  nucleons 
(neutrons  and  protons).  Let  us  Immediately  emphasise  several  Important  differences  between 
the  N-body  Hamiltonians  arising  In  Atosdc  and  Nuclear  Physics; 

1)  translation  Invariance  of  the  center  of  mass  In  Nuclear  Physics  (and  no  1-Body  terms) 

11)  very  different  2-Body  potentials  (In  Nuclear  Physics  the  potentials  have  very  short 
range ) 

111)  large  numbers  of  particles  (nucleons). 

* 

Service  P.T.N.,  Centre  d'Etudes  de  Bruybres-Le-Chatel,  B.P.  n‘’l2,  91680  Bruyeres- 
Le-Chatel  cedex. 

** 

Cereiaade,  Unlverslte  Parls-Dauphlne,  Place  de  Lattre  de  Tasslgny,  75775  Paris  Cedex  16. 
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We  will  come  back  on  these  differences  and  we  will  mention  others  such  as  the  use  of 
phenomenological  density-dependent  forces. 

From  the  the  mathematical  viewpoint,  these  differences  lead  to  equations  (Hartree-Fock 
equations,  HF  in  short)  of  a  completely  different  nature.  And  we  do  not  know  of  any 
reference  in  the  mathematical  (or  mathematic  physics)  literature  dealing  with  HF  equations 
in  Nuclear  Physics,  while  there  are  many  references  for  HF  equations  (or  at  least  Hartree 
equations)  in  Atomic  Physics  (see  for  example  E.  H.  Lieb  and  B.  Simon  [30],  E.  H.  Lieb 
[29],  P.  L.  Lions  [32]  and  the  references  therein). 

Our  goal  here  is  to  present  to  mathematicians  the  basics  of  Hartree-Fock  method 
(section  II  below)  together  with  the  known  mathematical  results  on  HF  equations  in  the 
context  of  Nuclear  Physics.  As  we  will  see,  many  problems  remain  by  large  open  and  we 
present  sometimes  model  simplified  problems  which,  hopefully,  should  preserve  the  same 
features  than  the  exact  HF  systems  of  equations. 

We  first  describe  the  HF  method  (section  II)  which  approximates  a  linear  problem  with 
a  single  unknown  function  in  large  dimensions  by  a  nonlinear  one  in  3  dimensions  with  a 
large  number  of  unknown  functions  (the  computational  advantage  being  obvious).  If  one  is 
interested  in  the  ground  state  of  a  nucleus,  the  resulting  problem  by  HF  method  is  roughly 
speaking  a  semilinear  vector  valued  minimization  problem  with  constraints  on  which  is 

translation  invariant.  This  is  typical  of  problems  which  can  be  analyzed  by  the  so-called 
concentration-compactness  method  (cf.  P.  L.  Lions  [33],  [34]).  We  explain  in  sections 
III  -  V  the  existence  results  we  can  obtain,  adopting  a  layered  presentation  to  cover  more 
and  more  realistic  problems  (from  the  physics  viewpoint).  However,  we  do  not  consider  in 
these  sections  the  possibility  of  spin-dependence  and  spin-orbit  forces  until  section  VI. 

In  section  VII  we  go  back  to  the  original  Hamiltonians  and  we  discuss  the  various 
approximations  including  Thomas-Fermi  classical  approximation.  Section  VIII  is  devoted  to 
the  search  of  solutions  of  HF  equations  with  symmetries  while  section  IX  is  a  very  small 
contribution  to  the  understanding  of  symmetry  breakings  of  the  nucleus.  In  section  X,  we 
describe  the  external  field  method  which  is  an  important  tool  for  the  numerical  computation 
of  the  ground  state. 
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Section  XI  is  devoted  to  various  considerations  on  time- dependent  Hartree-Fock  (TDHF 


In  short)  equations  such  as  the  orbital  stability  of  the  minima  of  HF  problems  and  the 
study  of  other  periodic  solutions.  From  the  mathematical  viewpoint  TDHF  equations  are 
systems  of  semlllnear  Schrfidlnger  equations. 

Finally,  the  last  section  (XII)  concerns  another  approximation  (somewhat  related  to 
the  HF  method)  known  as  the  Hartree-Fock-Bobolyubov  method  (HFB  In  short)  and  we  refer  to 
J.  Dechargi  and  D.  Gogny  [14],  J.  G.  Valatln  [47]  for  the  Physics  background  of  this 
method. 

There  are  Important  questions  related  to  HF  equations  that  we  will  not  consider  here 
neunely  the  question  of  numerical  analysis  of  HP  equations,  the  RPA  system  and  questions 
related  to  VIKB  approximations  when  goes  to  0.  We  hope  to  come  back  on  these  questions 
In  future  publications. 

Let  us  finally  mention  that  we  will  not  assume  any  knowledge  of  Quantum  Physics  from 
the  reader,  but  that  we  will  try  as  much  as  possible  to  keep  present  the  Physics 
motivations.  The  authors  would  like  to  thank  M.  R.  Dautray  for  bringing  HF  theory  In 
Nuclear  Physics  to  the  attention  of  the  second  author,  and  for  stimulating  their 
interdisciplinary  colloljoration. 
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II.  Presentation  of  H»rtree-Fock  method 

The  basic  object  we  consider  Is  a  A-body  Hamiltonian  that  we  denote  by  H,  where  A 
is  a  positive  integer:  in  Nuclear  Physics «  A  *■  N-t-Z  with  N  number  of  neutrons  and  Z 
number  of  protons.  The  precise  quantum  system  of  A  interacting  nucleons  is  supposed  to 
be  described  by  the  Hamiltonian  H 

(1)  j,  "i  " 

where  V  is  a  given  potential  (function  on  ft  is  the  Planck  constant,  m  is  the 

mass  of  the  nucleon  (we  neglect  here  as  usual  the  little  differences  of  mass  between 

m2 

neutrons  and  protons)  so  may  be  thought  of  as  a  given  positive  constant.  The  points 

x^( 1<i<A)  are  generic  points  of  and  the  notation  means  the  Laplacian  with 

respect  to  the  group  of  variables. 

The  Hamiltonian  H  is,  at  least  formally,  a  self-adjoint  operator  acting  on  the 
closed  subspace  of  L^((R^)*)  consisting  of  antisymmetric  functions  ♦  of  x  « 


,Xj^)  e  i.e. 


(2)  ^^^a(  1 ) '  *  *  * '^a(A)  ^  *  (**) )  ^ ^ ^ ^( x^ , . ,  •  ,x^) 

for  all  x^  £  (1<i<A)  and  for  all  permutations  a  of  {1,...,A}  where  |o|  denotes 

the  signature  of  o.  We  denote  by  H  this  subspace.  This  important  constraint  (2) 
corresponds  to  the  fundamental  Pauli  principle  and  is  due  to  the  fact  that  nucleons  are 
fermions. 

Before  going  further  in  the  description  of  HF  method,  let  us  point  out  that,  for 
physically  correct  Hamiltonians  H,  4  should  depend  on  spin  and  H  should  Incorporate 
spin-orbit  terms  and  density-dependent  3-body  terms.  We  deliberately  ignore  those  terms  in 
this  section  to  keep  the  ideas  clear  even  if  in  next  section  the  density-dependent  term  is 
incorporated.  Finally,  we  made  no  distinction  between  nucleons. 

Of  course,  the  interaction  is  mainly  described  by  the  choice  of  the  potential  V:  let 


us  mention  some  typical  examples  in  Nuclear  Physics 

(3)  V(x)  =  t  0e“'’l’'l  , 


t  Be"'’ I*  I  ,a,B€R,  u,v>0 


(4) 


V(x)  -  ae"''l’‘l(1/lx|)  +  e*"''l*l(1/|x|) 


0,6  «  R,  u,v  >  0 


or 

<5)  V(x)  “  o6j|  -  BASq,  o  e  r,  B  >  0 

Ml  theao  choices  (and  there  are  many  other a)  reepect  the  fundamental  character  of  strong 
interaction:  short  range  and  Intense  at  short  distances,  Notice  that  V  Is  spherically 
aymnatrlc  and  this  Is  also  a  general  feature  of  the  potentials  V  used  In  Nuclear  Physics. 
Observe  also  that  H,  H  are  Invariant  by  translation  of  the  canter  of  mass 


(6)  If  9  e  H,  Tjj*«*(x^+h,  Xjth, . . .  ,Xj^+h)  e  H  for  all  h  t  R^ 

(7)  “  ***'^h*>  ‘ 

Of  course,  one  wants  to  know  the  spectrum  of  H  and  Its  eigenfunctions.  In 
particular,  a  fundamental  role  Is  played  by  the  bottom  of  the  spectrum  which  Is  obviously 
given  by 

(8)  B  -  Inf{(H4,*)^2  /  ♦  «  H.  ‘ 

This  Is  the  so-called  ground  state  energy.  We  will  write  sometimes  to  recall  the 

dependence  on  the  number  A.  Of  course,  the  above  notation  Is  formal  since  (H4,$)  ^  la 

L 

not  defined  In  general  on  H  but  on  a  subspace  which  description  depends  on  V;  we  will 
Ignore  those  technicalities  In  this  section,  let  us  finally  mention  that  any  minimum  of 
(8)  Is  called  a  ground  state  (In  fact,  we  have  made  here  so  many  simplifications  that  one 
can  prove  there  exists  no  minimum  of  (8)  because  of  the  translation  Invariance  -  see 
section  VII). 

In  Nuclear  Physics,  one  has  to  deal  with  nuclell  for  which  the  number  \  of  nucleons 

Is  large  (up  to  240)  and  this  Is  why  the  direct  computation  of  (8)  is  hopeless.  Some 

approximation  Is  needed.  The  original  Idea  of  D.  Hartree  (  ]  was  to  consider  more 

A 

functions  ♦  (l.e.  tost  functions  ♦)  of  the  form:  4(x.,...,x,)  -  IT  9j^(xj^).  But 
clearly  this  choice  contradicts  the  antisymmetry  requirement  (2).  This  led  V.  Fock  [17 j 
and  J.  c.  Slater  [45]  to  a  better  choice  of  t  namely 
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(9)  «(x^,...,x^)  »— r")(-1)i  '  «  9a(i)(Xl)  -  det(»^(xj)) 

/ai  ®  1-1 

where  e^,...,ej^  are  A  functions  on  and  the  sum  Is  over  all  permutations  of 

{1,...,a}.  Such  a  choice  of  *  Is  called  a  Slater  determinant. 


Next,  to  check  the  normalization  constraint  of  *  In  (B),  we  see  that  It  Is  enough  to 


Impose 


Indeed,  we  then  have 


/  3  Si^iPjdx  -  for  1  <  1,  j  < 


3aI*I  ^  “  A1  ^  ^  ^  ^  3  ''c(l)*’‘l’’o'(l)^*l’'**l 

R  a,<j‘  1  R 


h-  I  i-’i 


a(l)o’(l) 


I  -  1  . 

a»a* 

Therefore,  the  HF  approximation  consists  in  replacing  E  in  (8)  by 

(11)  E„p  -  Inf{(H«,*)^2  /  *  ”  detCe^Cxj)),  /^j  -  5^^^} 

Observe  that  we  have  of  course 


Next,  It  Is  possible  to  rewrite  (H*,*)  3  quite  simply  when  ♦  Is  given  by  a  Slater 

£1 

determinant.  Indeed,  we  have  for  all  i  <  j 


R  a, a’ 


^  1  ^  3  ’a(k)‘’‘k’»^(k)‘’'k>’*’'k) 
k/l  R 


=  IT  - 

O  R 

-  X  5  /  3 1’*!  (*11^-1* 


'•  .  . 
*»'  V**  •- 


e-e^  -  -  -  >  '•••'  a***e*.»  * 


while 


~^l  !  |9^j^j(y)|2dxdy 

°  ir^tr 

R^xR^ 

And  we  obtain  the  following  expressions 

(13)  Ejjp  =  Inf  {E((p^  ,. . .  ,ip^)  /  (f^  £  L^(R^>,  /  ^  9^».dx  »  5^  for  1  <  1,  j  <  A} 


tA 


E(q  ,...,<p  ]  ^  I  /  <3x  + 

1=1  R^  ^ 


*  I  If  Ie.(x)|^  V(x-y) le  (y) l^dxdy  + 
T  1  *  3 


i<j  -3  3 

R  xR 


or  equivalently 


I  /  /  e,  (x)»  (x)  V(x-y)9  (y)q>  (y)dxdy 

i  ^  ^  3  *  J 


R^R^ 


(  15  )  Et®^,  . .  .  =  -jjj-  y  J  i  I  ^dx  +  -i  T  /  /  I  *^(x)  I  ^V( x-y)  |  »  ,  (y)  |  ^dxdy 


R^R^ 


y  y  /  /  ».(x)e  (x)  v(x-y)e  (y)»  (y)dxdy  . 


r’xr^ 


the  second  term  is  often  called  the  direct  term  while  the  third  one  is  called  the  exchange 

A 

term.  We  will  also  often  denote  by  t(x)  the  density  of  Kinetic  energy  T  *  \  |V9^(x)| 

A  ^  i-1 

P(x)  the  density  P  *  I  '(p*(x)’^  and  p(x,y)  the  density  matrix  p(x,y)  *» 

A  i-1 

r  * 

®.(x)®j(y)«  Observe  that  we  may  simply  write  F  as 


E[<P,  » •  •  • »».  ]  “  T"  ^  ^  !  P(x)V(x-y)p(y)dxdy  + 

1  A  2lTt  '3  ^  ^  ^ 


(16) 


R^xR^ 


“7  /  /  V(x-y) |p{x,y) l^dxdy  . 

R^R^ 


The  HF  minimization  problem  (13)  Is  nonlinear,  nonconvex  In  general,  with  constraints 
and  Is  invariant  by  translations  (translating  at  the  same  time  all  ,  rotations  In 
R^  and  by  unitary  transforms  of  (e^,...,^^^)  in  C^.  In  fact.  In  the  examples  given  In 
next  sections,  the  HF  minimization  problems  will  l9e  slightly  different  (but  still  with  the 
same  general  features);  an  additional  nonlinear  term  will  be  Included  In  E  corresponding 
to  a  density-dependent  nonlinear  3-body  term  In  H,  We  will  come  back  on  the  realistic 
H  being  used  in  Nuclear  Physics  in  section  VII  while  we  analyze  in  section  IV  -  V  HF 
problems  like  (13)  deduced  from  various  examples  of  these  realistic  Hamiltonians  by  the 
method  described  above. 

The  Euler- La grange  equations  corresponding  to  the  minimization  problem  (13)  may  be 
written  as 


h 

2m 


+  (0*V)9^  -  /  j  g^(y)  [-j  V(x-y)  +  7  V(y-x)  ]  p(x,y)dy 


R 

1  e 

j 


in  R 


for  some  hermitian  matrix  Lagrange  multipliers.  Now,  observe  that  if  U  is 

unitary  and  diagonalizes  then  (e^,...,9^)  “  U(»^ , . . . ,  Sj^)  is  still  a  minimum  if 

(9,,...9  )  minimizes  (13).  And  (».,••.,*.)  now  solves 
'  «  1  A 

(  17) 


-  ^  A9i  +  (0*V)(p^  -  =•  e^9^  in  R^ 

for  some  constant  ej^,  where  K  Is  the  operator  defined  by 

K(p(x)  =  /  j  9(y)(7V(x-y)  +  7  V(y-x)  ]  D(x,y)  dy  . 


In  particular,  the  constants  e.,...,e.  are  eigenvalues  of  the  operator 


(-  A  (P*V)  -  K}«  The  syetem  of  equations  (16)  Is  called  HF  equations*  We  also  remark 

2 

that  are  eigenvalues  of  the  above  operator  H  “  -  A  +  (p*V)  -  K. 

Notice  also  that,  at  least  in  all  examples  considered  below,  these  eigenvalues 
are  non  positive. 

We  conclude  this  section  by  a  brief  discussion  of  the  validity  of  HF  method;  notice 
that  it  is  an  approximation  of  the  *true“  problem  (8)  and  that,  a  priori,  it  gives  only  a 
bound  from  above  of  the  ground  state  energy  E  (recall  (12)).  On  the  other  hand,  there 
are  various  reasons  to  use  it  and  thus  study  it:  first  of  all,  it  gives  good  numerical 
results  and  there  are  almost  no  substitutes  to  compute  the  ground  state  on  E.  A  more 
"scientific"  reason  is  its  asymptotic  validity  as  A  ■*•  +«  (for  general  V)  as  proved  by 
E.  H.  Lieb  and  B,  Simon  [30],  [31]:  we  will  come  back  on  this  point  in  section  VII. 


III.  A  model  case 


I 


To  give  an  idea  of  the  type  of  HF  problems  idilch  are  encountered  in  Nuclear  Physics / 
we  will  consider  In  this  section  a  vary  simplified  problem:  we  build  a  scalar  problem 
(A  >  1)  with  Skyrme's  Interaction  as  In  D.  Vautherln  and  D.  M.  Brink  [48]  skipping  the 
spin  dependence.  In  the  next  section,  we  will  consider  the  general  case  of  HF  problesw 
with  Skyrme's  Interaction  but  without  spin  dependence.  The  functional  (16)  becomes  in  this 


■  to  ^  -  I  /  jlel^dx  +  I  ;  jlel^dx  + 


+  4/  / 


where  a,  8,  y,  S  are  constants  such  that:  a  >  0,  8  >  0,  y  >  0,  y  +  fi  >  0.  This 

functional  corresponds  (essentially)  to  a  potential  v  of  the  form  (5)  with  the  additional 
3  f  t  1 6 

term  (j  /  <p  dx) . 

Then  the  HF  minimization  problem  becomes 
(19)  I  •  inf{E(»)  /  e  X,  /^lepdx  -  1} 

where  the  minimizing  class  X  Is  defined  by 

X  -  {9  «  H^R^)(*),  /  Jsl^lVel^dx  <  ">)  . 

ST 

Before  stating  our  main  existence  result,  let  us  recall  that  a  minimizing  sequence  (p^) 

Is  a  sequence  (9„)  In  X  satisfying 

’  '  '‘V  5  ^  • 

We  have  the 

Theorem  III .  1 ;  For  every  minimizing  sequence  (Tn^n  minimization  problem  (19)  one 

can  find  In  R^  such  that  the  new  minimizing  sequence  relatively 

compact  in  H^R^)  If  and  only  If  I  <  0.  In  particular,  if  I  <  0  there  exists  a 


(*)  h’(R^)  =  [,  :  L^(R^),  L^(R^)} 


.  'V'  A  *  •_ 


.  *  .  •  m  s"'  a**  •  •  **'_r'*  •'*  •  •  *  * 


'  .  ,  •  .'  a'  H  - 

•-V  \  ^  S  ^  * 


•  ■  •  ■  I .  V  .■ '!•  'j ■^"j.  vu ",  'A'l'  'a v»" w »"» 


l» 


:-v>.'--. 


(20) 


-  Inf{E(<»)  /  »  e  X,  /  jUl^dx  -  X) 


where  X  >  0.  Then,  applvintj  the  arguments  of  (33l  ,  we  deduce  that  any  minimizing  sequence 
of  (19)  is  relatively  compact  up  to  a  translation  if  and  only  if  I  the  following 
condition  holds 

(S.1)  I,  <  !„  +  I,.„  ,  Va  €  (0,1)  . 
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minimum  of  (19).  In  addition,  I  <  0  if  and  only  if  a  >  Og  where  Og  is  a  positive 

^2  k2 

constant  depending  on  B,  y,  6,  which  goes  to  0  as  goes  to  0.  If  a  <  Jr  ■  ir«^ 

there  are  minimizing  sequences  converging  to  0  in  I^(B^)  for  2  <  p  <  ».  Finally  if 
a  <  Og  there  is  no  minimum  of  (19). 

Remarks ;  i)  Further  properties  of  minima  of  (19)  are  given  below. 

ii)  Scalar  problems  lilce  (19)  have  already  been  studied  by  several  authors:  we  will  only 
mention  the  trorlis  by  W.  Strauss  [46],  C.  V.  Coffman  [13],  H.  Berestycici  and  P.  L.  Lions 
[5],  P.  L.  Lions  [40].  The  methods  in  these  wor)ca  yield  only  the  existence  of  a  minimum 
if  I  <  0,  using  a  symmetrization  argument  which  is  outlined  in  the  proof  of  Proposition 
III. 2  below  and  which  no  longer  applies  to  more  realistic  problesw  such  as  the  ones  studied 
in  the  following  sections. 

iii)  In  fact  it  is  possible  to  treat  the  case  ft  •  0:  in  that  case  every  minimizing 

sequence  in  the  class  (e  e  L^(R^);  (He  e)^,  (Im  9)^  e  H^(R^);  f  ls)^|7e|^dx  <  ■}  is 
relatively  compact  say  in  L^  up  to  a  translation  and  there  exists  a  minimum  which  is  the 
limit  of  the  minima  of  (19)  as  fi  0.  □ 

Before  proving  Theorem  III.I,  we  prove  the 

Proposition  III. 2:  Assume  that  I  <  0.  Then  there  exists  a  minimum  of  (19)  which  is 
spherically  symmetric,  positive  and  decreasing  with  respect  to  r  •  |x|. 

Proof;  In  view  of  Theorem  III.I,  there  exists  a  minimum  p  of  (19),  Then  considering  the 
spherical  nonincreasing  rearrangement  of  {vl,  one  chec)cs  easily  that  E  is  decreased  and 
thus  a  minimum  with  the  above  properties  is  found.  □ 

We  now  turn  to  the 

Proof  of  Theorem  III.I;  We  are  going  to  apply  the  concentration-compactness  arguments  (cf. 

P.  L.  Lions  [33] ).  Hence  we  introduce 


i. 


•V- 


•  V, 


We  will  not  redo  the  proof  In  [33]  but  we  will  only  make  a  few  formal  obaervatlons  in 

order  to  explain  the  role  of  The  main  difficulty  in  the  above  statement  is  the  "if 

part:  in  [33],  it  was  proved  that  if  a  minimizing  sequence  is  not  relatively  compact  up  to 

a  translation  then  roughly  speaking  it  breaks  at  least  into  two  parts  which  are  essentially 

supported  in  two  disjoint  closed  sets  whose  distance  goes  to  Let  us  denote  those  two 

parts  by  so  e  ~  we  may  assume  that  /  ,(e')^dx  ♦  a, 

nn  n-^n  n  n  n 

r  2  2  * 

I  (e  )  dx>  1-a.  The  above  dichotomy  then  yields 
n  n 

R  12 

I  -  lim  E[e  ]  >  lim  £[«  ]  +  lim  E[e  ]  >  I  +  I,  ^ 
n  '  —  n  "  '  n  a  1-a 
n  n  n 

and  if  (S. 1)  holds  this  is  not  possible. 

The  arguments  of  [33]  apply:  the  only  modification  consists  in  checking  that  if 

is  bounded  in  H^(R^),  /  3Wnl  °  (indep.  of  n)  and  converges  weakly 

R 

in  to  some  e  then 

n  R  R 

R  R 


The  proof  of  this  claim  is  a  simple  consequence  of  Lenraa  III. 3  which  is  stated  and  proved 
after  the  proof  of  Theorem  III.1. 

We  next  show  that  (S.1)  is  equivalent  to  I  <  0  and  that  one  has  always  <  0: 

Indeed,  let  9  e  P(R^)  such  that  /  ^lel^dx  =  X  and  let  9g(x)  •  o  for  a  >  0. 

R 

Obviously 


OR  OR 

+  \  {|-  /  +  "T  tj"  /  ,ll>|^|ViI>l^dx  +  ^  /  jlvlel^l^dx} 

R 


6  '6  ■'  3' 

O  R 


5  ‘4  ^  3' 

O  R 


and  thus  letting  o  +  +■,  we  prove  that  <  0.  Next,  by  a  similar  scaling  argument  we 


see  that  for  all  9  >  1,  X  >  0 


^0X  ”  ^  4  /  ^  I  jl'fl^d-'tl  + 

R  R  R 

+  /  jl»l^|7»|^dx  +  ^  /  jlvl»|^|^dx)/*£X,  /  jUl^dx  -  X}  . 


Now  we  claim  that  If  <  0  the  Infimum  in  ly  may  be  restricted  to  those  s  satisfying 
2 

Ats)  ”  ^  j  |7s|  ^dx  +■  /  j  |s|  ^|7s  I  ^dx  +  /  j  l7  j  f  |  ^)  ^dx  >  v  for  some  u  >  0.  Indeed  if 

R  R  R 

there  were  a  minimizing  sequence  of  Ij^  that  ACSj^)  *  0,  then  by  Sobolev  embeddings 

9^  ♦  0  in  IiP(R^)  for  2  <  p  <  6  (in  fact  p  <  12  here)  and  so  /Is^l^dx  ♦  0.  But  this 

means  that  >  0  and  our  claim  is  proved.  Hence,  we  deduce 

Igj^  <  8  lnf{E(9l/9  e  X,  A(9)  >  V,  /  ^Isl^dx  -  X}  -  81^^ 

R 

and  this  inequality  holds  if  I;^  <  0.  Then,  a  straightforward  argument  proves  that  (S.1) 
is  equivalent  to  I  <  0. 

Observe  also  that  the  above  scaling  argument  shows  that  if  I  -  0,  there  is  a 
minimizing  sequence  (namely  Sg  as  a  *  +•)  which  converges  to  0  in  for  2  <  p  <  •. 

He  next  discuss  the  inequality  I  <  0.  It  is  obvious  that  I  <  0  for  a  large 
enough.  So  let  us  denote  by  “  ®o  least  positive  constant  such  that 

I  <  0  for  a  >  Og  (observe  that  I  is  nonincreasing  with  respect  to  a).  He  have  to 

prove  that  Og  >  0  or  in  other  words  that  I  >  0  for  a  small  enough.  But  using  Sot>olev 

and  Holder  inequalities  we  find  for  geX,  / 

R^ 


(/  jlsl^dx)’''^  ♦  |(/  jlsl^dx)  . 

R  R  R 


A  simple  study  of  the  function  of  one  real  variable  Cg  ^  t'^^  -  t'^^  6^  ^  proves 

that  I  •  0  for  a  small  enough. 

The  fact  that  Og  goes  to  0  as  — ♦  0  can  be  seen  from  the  expression  of  E[9„l 

V  Ttt  o 

given  above. 

Finally  let  a  <  a  ,  if  there  exists  a  minimum  of  I  for  a  we  can  test  the 


m'’  m  •  •  w*  •  •*'  •  **•  4 


V.".-. • 

•  •  .%  -w'-  /  . 


minimization  problam  (19)  for  a  <  a  <  with  this  minimum  and  this  glvas  a  nsgatlvs 
value  for  I  contradicting  the  definitl''n  ''  Oq.  a 

I^Hma  lll.3»  Let  ^  0  converge  In  L^  to  eome  f  >  0,  let  g,^  converge  wea)cly  in 
L^  to  some  g  and  assume  that  for  19^1^  bounded  In  L^.  Then  f|g|^  e  and  we 


have 


ii®  /  >  /  f|gl‘ 

n 


Proof;  Let  M,  R  e  (O.o).  We  Introduce  a  symmetric  convex  function  satisfying 

0  <  ^j^(z)  <  Ul^>  ~  for  |z|  <  R,  is  Lipschltz.  It  is  clearly  enough  to 

prove  that  we  have 

li»/*nl9nl^*  /  «  • 


n  (f<M) 

In  order  to  do  so  we  first  observe  that,  without  loss  of  generality,  we  suty  assoae  that 


2 

(g„)  converges  weakly  in  L  to  some  h  which  satisfies  h  > 


Next,  we  remark  that  we  can  conclude  if  we  prove  that 


/  (f-f  >  ♦  0 

(f<M)  »  R  “  » 


But  this  integral  is  easily  bounded  for  all  6  >  0  by 

« f  >«)VV  * 


<  «  /  |gj2  +  C(M,R)  ,5,|gJ 


<  C5  +  C<MfR){me«8  (f-f  >  6)}'^ 


<  C«  +  C(M,R)«”  f  I  f-f 


and  we  may  conclude 


IV.  8]cyrma*«  lnt«raetlon  without  spin 

In  all  this  ssctlon  we  will  consider  only  the  so-called  Skyree's  Interaction 
thus  following  the  approach  by  D.  Vautherln  and  D.  H.  Brink  [48].  To  slapllfy  the 
presentation  we  skip  the  spin  dependence  and  we  refer  to  section  VX  for  the  coo^lete 
problems. 

We  begin  this  section  by  a  special  case  which  corresponds  to  the  slsqpllfled  situation 
where  no  differences  are  made  between  neutrons  and  protons  (N>Z,  no  Coulomb  Interaction 
between  protons).  In  fact,  even  If  Isospln  If  "fixed*  for  neutrons  and  protons  at  the 

level  of  the  original  X-body  problem,  one  can  allow  In  the  Slater  determinant  wave- 

functions  which  depend  on  the  Isospln  and  then  one  Is  also  led  to  probleu  of  the  following 


form.  In  those  cases  the  functional  (16)  becomes 


(21) 


®t^1 . '  f  *  I 


A  ^2 

where  t  -  J  1  »  “  I  |■■,l  •  The  constants  o,  6,  y,  S  satisfy 

1-1  1-1 

(22)  a  >  0,  8  >  0,  «  >  0,  0+Y  >  0  . 

As  In  the  preceding  section,  only  the  last  term  T  /  ,  Is  not  an  obvious  consequence 

* 

of  the  HF  method  as  described  In  section  II:  Indeed  this  tsrm  (and  analogous  terms  In  this 
section  and  in  the  next  one)  comes  from  a  3-body  term  which  Is  equivalent  to  a  2-body 
density  dependent  term  In  the  Hamiltonian  H.  We  will  come  back  on  this  point  In  Section 
VII. 

The  HF  minimization  problem  may  then  be  written  as 


I  -  Inf{E(e^,...,e^]  /  €  h’(R^),  /  j  pTdx  <  », 


(23) 


/^3  v,f;dx 


6^^  for  1  <  1,J  <  A}  . 


Before  stating  our  main  result  on  (23),  we  need  to  Introduce  a  few  notations:  let 
M  -  (m^^j)  be  a  nonnegative  hermitlan  matrix,  we  Introduce  the  following  ndnlmleatlon 


problem 


Ijj  -  /  ?^  e  /  3  PTdx  < 


(24) 


3  ’I’j*** 


for  1<i,J<A} 


Then,  we  observe  that  since  B  Is  Invariant  under  unitary  transfonts  of 
for  any  unitary  matrix  U  in 


then 


(25) 


I  -  I 
H  0  ’MO 


so  we  may  choose  0  so  that  0~^NU  Is  diagonal  and  If  (m^,...,mj^)  are  the  eigenvalues 
of  M,  Mg  »  dlag(m^, . . .  ,m,^)  Is  the  diagonal  matrix  with  (m^,..,,mj^)  as  diagonal  entries 
then  Ijj  ”  Ij(^  and  we  will  denote  Ijj  *  »  I(m^, . .. ,mj^).  With  these  notations,  I 

*  1(1, ...,1)  where  1-  (^ijJ* 

Finally,  «re  will  say  that  a  sequence  ...,  in  the  minimizing  set  is  relatively 

compact  up  to  a  translation  If  there  exists  y"  in  f?  such  that  (i^,  •••,e||^)  ■ 

(»"(*+y"), . . . ,f^( •+y") )  -  which  la  still  a  minimizing  sequence  -  is  relatively  compact 

In  and  |V  p"]^,  p"  r"  are  compact  in  (with  obvious  notations). 

Theorem  IV. 1 i  1)  The  Inflmum  I  e  (-“,0]  and  for  all  R  <  •  there  exists  <  •  such 

that  /  3  p  +  T  +  pT  dx  <  Cjj  if  K(f,,...,f,^l  <  R,  /  3  “  ®lj*'^l,j'' 

/  3  PT  dx  <  •. 

« 

11)  Every  minimizing  sequence  of  the  problem  (24)  Is  relatively  compact  up  to  a 
translation  If  and  only  If  the  following  condition  holds 

(26)  1  <  Kre^ , . . .  ,mj^)  +  1(  , . . . ,1*mj^)  for  all  (re^,...,m^) 

such  that  0  <  m.  <  1  for  1  <  1  <  A,  0  <  ^  m.  <  A.  Of  course,  if  (26)  holds  then  there 

1 

la  a  minimum  of  ( 23 ) . 

>.2 

111)  If  a  <  Ug  where  Ug  Is  some  positive  constant  depending  on  A,  ,  B,  o,  $>  then 
I  >  0  and  there  Is  no  minimum  of  ( 23 ) . 

Remarks;  1)  In  general,  we  do  not  know  how  to  check  (26).  The  answer  seems  to  be  highly 


dependent  on  A  In  view  of  the  numerical  computations  which  have  been  performed.  In  any 
case,  checking  conditions  (26)  when  numerical  computations  of  ground  states  appears  to  be  a 
good  test  since  (26)  means  a  certain  stability  of  the  absolute  minimum. 
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11)  In  fact,  as  seen  below  from  the  proof  which  again  relies  on  the  concentration- 
compactness  arguments  [33],  the  concentration-compactness  method  not  only  shows  the 
necessity  and  sufficiency  of  (26)  but  also  predicts  what  can  happen  on  minimizing 
sequences.  Imt  >.3  give  a  few  examples)  1)  suppose  I  -  0,  then  there  are  minimizing 
sequences  converging  to  0  In  for  p  >  2  and  the  density  vanishes  (In  the  sense  of 

[33]  ), 

2)  suppose  (to  simplify)  that  there  exists  a  unique  set  of  values  (m^,...,m^)  such 
that  0  <  mj^  <  1  for  1  <  1  <  A,  0  <  I  m^^  <  X,  I  -  I(m^ , . . .  ,mj^)  +  I(  1-m^ ,  • . . ,  l-m^^)  while 
(26)  holds  for  all  (m^,...,mj^)  ^  (m^ , . . .  ,mj^) .  In  fact,  this  over  simplification  Implies 
m^>.  ..«m^  but  we  will  Ignore  thle  for  the  sa)ce  of  the  argument.  Two  cases  may  occur:  the 
slmpleat  one  la  when  the  two  minimization  problems  I(m^ ,. . .  ,m^),  I(  1-m^ , . . . ,  l-m^^)  satisfy 
the  analogous  of  the  aubaddltlvlty  conditions  (26).  Then,  there  are  minimizing  sequences 
of  (23)  which  are  not  relatively  compact  up  to  a  translation  and  any  such  sequence 
(g!|,...,e^)  brea]ca  Into  two  parts: 

e“  -  ♦"  +  xj  »  1  <  1  <  A 

where  <)",  x"  Are  relatively  compact  up  to  a  translation  and  are  minimizing  sequences  of 

I(m^ ,..  .,i^),  I(  , . . . ,  and  thus  (extracting  subsequences  If  necessary)  converging 

to  mlnlM  of  these  problems.  In  addition,  roughly  spealclng  the  distance  between  the 

supports  of  J  and  J  lx?|^  goes  to  “  as  n  ♦  •*.  The  second  case  concerns  the 

1  _  _  1  ^  _  _ 

situation  when  I(m, ,...,m. )  (or  I ( 1-m. , . . . , 1-m,) )  does  not  satisfy  the  analogue  of 
1  A  I  A 

(26):  then  tre  may  continue  the  above  argument  and  In  turn  i|>"  can  break  Into  two 
pieces.  If  we  knew  completely  the  function  I(m^ , . . . ,mj^) ,  It  would  be  possible  to 
determine  completely  the  behavior  of  islnlmlzlng  sequences:  vanishing,  dichotomy  Into  n 
parts  converging  to  minima  of  subproblenw ,  dichotomy  Into  n  parts  with  (n-l)  pieces 
converging  to  minima  of  subproblems  and  one  piece  vanishing. 

Ill)  Again,  we  can  treat  as  well  the  case  ))  •  0  and  the  analogue  of  1)  holds. 

a 
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At  this  stage.  It  Is  Interesting  to  briefly  explain  how  all  the  above  phenonena  are 
related  to  various  physical  situations  which  essentially  depend  on  the  set  (N,Z) 
characterizing  to  numbers  of  neutrons  and  protons.  In  fact.  It  Is  experimentally  observed 
that  the  existence  of  a  nucleus  crucially  depends  on  the  set  (N,Z)  as  we  explain  now  with 
the  help  of  Figure  1  below  taken  from  M.  Effer  [IS].  In  the  plane  ()I,Z)  the  stable 
nuclei  (Infinite  llfe-tlme)  are  Indicated  In  black.  Admitting  that  the  A-body  Hamiltonian 
and  the  Hartree-Fock  approximation  correctly  represent  the  reality,  these  values  of  (M,Z) 
would  correspond  to  "nice”  minima  In  the  problems  we  are  considering  here  and  below  (the 
strict  subaddltlvlty  Inequalities  should  hold  for  such  values).  The  dotted  grey  zone 
corresponds  to  unstable  nuclei  which  are  known  today  and  whose  llfe-tlme  may  vary  In 

r 

between  10^^  years  and  some  mllll-seconds.  Let  us  mention  that  a  little  more  than  2000 
nuclei  are  known;  about  300  exist  In  nature  while  1900  were  "built”.  Bet%raen  2000  and 
4000  more  nuclei  are  expected  to  exist  (mostly  unstable).  Finally,  the  majority  of  nuclei 
currently  observed  In  nature  (263  out  of  287)  are  stable.  The  white  zone,  delimited  by  two 
lines,  correspond  to  (unstable)  nuclei  which  are  to  be  discovered.  For  HF  problems,  those 
unstable  nuclei  correspond  to  minimization  problems  where  the  strict  sub-addltlvlty 
Inequalities  do  not  hold  and  minimizing  sequences  break  Into  severl  "compact"  pieces  (see  a 
precise  example  below).  The  two  lines,  the  so-called  "drlp-llnes" ,  beyond  which  no  nuclei 
are  expected  to  exist,  are  precisely  associated  with  the  loss  of  exactly  one  neutron 
(Sjj  =  0)  or  one  proton  (Sp  =  0).  In  our  context,  this  situation  would  correspond  to  the 
case  when  (m,,...,m.)  -  (1,0,. ..,0)  in  Remark  11)  above  l.e.  minimizing  sequences  break 
Into  two  parts:  one  which  Is  "compact"  and  converges  up  to  a  translation  to  the  minimum  of 
a  I(0,1,...,1)  while  the  other  part  vanishes.  The  zone  beyond  the  drip  lines  should 
correspond  to  similar  phenomena  where  minimizing  sequences  break  into  several  pieces  one  of 
which  vanishes. 

In  order  to  Illustrate  the  situation  concerning  the  unstable  nuclei  we  shall  restrict 

144 

ourselves  to  two  examples.  The  first  one  concerns  the  nucleus  Nd ( Z“60 , N“84 )  whose 

lifetime  la  quite  long  (about  2  «  10®  years).  This  nucleus  Is  unstable  and  eventually 
decays,  emitting  an  alpha  particle  (elementary  nucleus  composed  of  2  neutrons  and  2 
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protons),  into  ths  two  stable  subaysteas  and  a  and  one  writes 

144  140 

N.  ♦  X  +■  a  . 

60  58  • 

The  other  exaaple  is  provided  by  the  nucleus  Pu  i«hich  spontaneously  fl^e  ons  Into  two 

stable  nuclei 

240  134  106 

Pu  ♦  T  ♦  . 

94  52  •  42  ° 

In  our  context  (once  aore  admitting  the  models  are  good  enough  to  reproduce  theee  fissions 
and  this  seems  to  be  the  case  In  view  of  HP  numerical  computations)  this  obviously 
corresponds  to  minimizing  seguences  breaking  into  two  parts  which  converge  (up  to 
translations)  to  minima  of  appropriate  subproblems.  In  fact,  if  we  were  to  use  a  nuclear 
force  realistic  enough  and  If  we  knew  completely  the  functions  I(m^ ,  • . .  ,a|^)  we  would  be 
able  to  predict  all  unstablllty  patterns  (or  confirm  the  numerical  computations  at  least  ..)■ 
We  wish  to  conclude  these  physical  considerations  by  Indicating  that  HP  minimisation 
problems  (with  possibly  the  extension  to  HFB  problems  •  see  section  XII)  lead  to  numerical 
computations  which  reproduce  quite  well  at  least  parts  of  the  dlagr^1m  below  (stable  nuclei, 
soma  unstable  ones,  drip  lines...  )t  the  rmttrlction  being  essentially  due  to  the 
difficulty  of  solving  numerically  these  problems.  And  we  refer  to  J.  P.  Berger,  N.  Glrod 
and  D.  Gogny  (6),  M.  Glrod  and  B.  Graramatlcos  (211,  0.  Vautherln  and  D.  M.  Brink  [48],  P. 
Quentin  and  H.  Flocard  [44],  J.  Negele  [41)  and  (42)  (and  the  references  given  therein)  for 
various  extensive  computations,  knother  observation  consists  In  remarking  that  for 
unstable  nuclei  In  fact  several  different  f ragswntatlons  are  often  possible  with  one  being 
more  probable  (statistically)  and  these  various  choices  could  be  related  to  dichotomies  of 
minimizing  sequences  corresponding  to  values  of  m^,...,mj^  strictly  between  0  and  1. 
Finally,  we  wish  to  warn  the  Interested  reader  that  the  above  considerations  Indicate  that 
strict  subaddltlvlty  Inequalities  may  be  very  hard  to  chock  and  in  addition  should  depend 
In  a  sensitive  way  on  A  (or  (N,z)  for  problems  below...). 


V  ■  -  '  *  •  ►  '  •  \ 


Figure  1 


Proof  of  Theorem  IV.  It  it  Is  enough  to  prove  the  existence  of  Cj^.  Remark  that 
I  PT  +  ■^  |Vp1^  -  ^  PT  ♦  tJ  Re(e^7<f*)  pT  >  0 

2 

01  2  ^2  ^  S^Y 

and  j  p  <  5  P  +  C(a,«)p  so  E [s, , . . .  ,Sj^l  >  /  3  •j;’  f  ♦  •'^4  )Pf  "  C(a,6)  and  i)  is 

R 

easily  deduced.  Again  part  ii)  of  the  above  result  is  a  direct  application  of  the 
concentration-compactness  argument  ([  ]).  We  will  not  give  the  proof  but  instead  we  will 

explain  the  main  idea  used  to  prove  the  sufficiency  of  (26).  If  (s"t...>e]^)  is  a 
minimizing  sequence  of  (23),  then  we  apply  the  lemma  below  (proved  in  (331)  with  the 
prolMibility  Pjj  on  R^  whose  density  is  ^  p"  (i.e.  the  density  in  Nuclear  Physics 
terminology! ) 

Lemma  IV.  2;  Let  (P^),,  be  a  sequence  of  probability  measures  on  R«.  Then  there  exists 
subsequence  that  we  still  denote  by  P^^  such  that  one  of  the  following  properties  hold 


i)  (compactness  up  to  a  translation)  3y''  €  Ve  >  0,  3R  <  « 

Pn(B(y",R))  >  1-e 

ii)  (vanishing)  WR  <  Sup^  P„(B(y,R))  ♦  0 

yeR 

ill)  (dichotomy)  3o  e  (0,1),  VC  >  0,  VM  <  “,  3r^  >  M,  ^y^^  c  if* i  3Rjj  ♦  *  such  that 

l*’n'®<yn'’^o’>  -  “I  <  E,  lPn(B(y„,Rn)‘=)  -  (1-a)|  <  e.  o 

1  A 

If  Pjj  or  Pjj  vanishes  (case  ii))  then  (see  [33])  <p^,...,q>J{  converge  strongly  in 

I/P(R^)  to  0  for  2  <  p  <  6  (actually  <  12)  and  thus  I  =  lim  E[iPj|,  . .  .  ,<p^)  >  0.  Since 

n 

one  checits  easily  by  a  scaling  argument  as  In  the  proof  of  Theorem  III.I  that  I|,j  <  0, 

this  means  1^0  and  it  contradicts  ( 26 ) . 

If  dichotomy  occurs  (case  ili))  then  we  translate  (<p^,  . .  •  ,?„)  by  y"  and  roughly 

spealcinq  we  split  these  functions  into  their  "restrictions"  to  B(y",Rq)  and  to  B(y",Rj^)‘^ 

and  we  denote  by  (i(i^, .  . .  ,i|<^) ,  ()(^,«--,Xn)  the  two  parts.  We  may  assume  that 

f  dx  *  m.  j  dx  ♦  6 .  .  -  m,  ^ 

"  "  n  ij  r3  "  "  "  ^3  ij 

for  some  hermltian  matrix  such  that  (essentially)  ^  =  oA.  The  contradiction 

with  (26)  is  obtained  by  remaricing  that 

I  =  11m  E(s  \  . ..  ,<p^l  >  1^  E(\)i\  .. .  +  lim  E(x  \  . . .  ,X^1 

n  n  n 

»  "m  ^  • 

Therefore,  if  (26)  holds  then  automatically  we  are  in  case  i)  and  we  conclude  as  in  [33] 
provided  one  observes  that  since  1  is  not  assumed  to  be  positive  there  is  a  little 
difficulty  to  pass  to  the  limit  which  is  solved  by  the 

Lemma  IV. 2;  Let  converge  wealcly  in  H^(R^)  to  Assume  in  addition 

that  is  bounded  in  l’(R^).  Then  pt  €  L^(R^)  and 

(27)  lim  /  ,PTdx>/,PT  dx,  lim  /  .  P  T  -  4  i  "  T  l^P  I  . 

^3  n  n  „3  -  _3  n  n  4  n  3  4 

n  R  R  n  R  R 

Proof !  The  first  part  of  (27)  is  a  consequence  of  Lemma  III. 3.  The  second  part  will 
also  be  after  a  few  considerations.  We  introduce  the  nonnegative,  convex,  quadratic 
functional  for  all  z  e  C*,  0  »  (9^,...,Sj^) 


Q(z,V0)  =  — 


(l  +  lz^) 


(I  IVP,!^) - flT" 

i  >•  (1  +  lz|2)  i  >■  1 


■■■ 


Observe  that  ■  (1+P)  [px  -  —  7p 

4 

Now,  if  we  set  )^2  ,  is  bounded  in  L^(E^).  And 

we  may  assume  that  g^^  converges  weakly  in  L  to  some  g.  If  we  show  that 

V  1  A 

g  >  {<5(0, 7s)}  2  where  9  “  (s’,...,s"),  then  applying  Lemma  III. 3  with  f^  ” 

conclude  easily  showing  the  second  half  of  (27)  on  any  bounded  domain  of  (this  is 

1  2 

enough  since  pT  -  —  |7p|  >0). 

1/ 

Since  Q(2,70)  2  is  convex  for  all  z,  it  is  clearly  enough  to  show  that  for  any 


*  £  D^sr) 


!  ,  iI-Iq'^2  (0  (U,70)|dx  0  . 

J  n  n  n 


1/ 

Since  Q  2  is  continuous  on  bounded  sets,  we  introduce  the  local  modulus 


0)^(6)  =  sup{  Iq'^2  (2,p)  -  Q''2(z+h,p)|  /  IpI  <  R,  |z|  <  R,  |h|  <  6} 


J/2 


for  all  5  >  0,  R  <  “ 

and  we  split  the  above  integral  into  integrals  over  several  sets  that  we  bound  as  follows: 

•  on  >  R  or  1 70 1  >  R  or  l0^!  >  R  or  |0|  >  R),  the  integral  is  bounded 

by 

C  /,  iK|70  I  +  |70|)dx  <  C(|(|)  meas(-)'^2  < 

I  * )  n  K 

•  on  /  |l*l  ^  R>‘  l•*n”**l  ^  integral  is  bounded  by 

C  /  iK|70  |  +  |70|)dx  <  C(i(-)  raeas(|0  -0)|  >  6)^2  ^  q 

(|0„-0|>«) 

•  on  {1 70  I,  1 70],  |0  ],  |0|  <  R;  |0  -0|  <  5),  the  integral  is  bounded  by 

n  n  n 

a)j^({)  /  j  i(idx  =  C 

R 

This  enables  us  to  conclude  easily.  Another  possible  proof  (communicated  to  us  by 
H.  Breris)  is  the  following:  since  ”  J  |7pj^|^,  PT  -  J  |7pj^  are  nonnegative,  it  is 

enough  to  show  that  for  all  M 

lim  ,f  p  T  -  |7p  l^dx  >  /  pi  -  T  |7p|^dx 

n  n  (  |9^|<m)  n  (  U^|<m) 

i=1  i=1 


Now,  by  Eqrov's  theorem,  for  all  e  >  0  there  exists  a  set  E  such  that  its  complement 
has  measure  less  than  £  and  converges  uniformly  to  on  E  for  all  i.  Denoting 


by  F  »  E  n  (  n  (|9  I  <  M) ) ,  It  is  enough  to  show  that 
i-1 


Itm  /_  P  T  -  x  |Vp  l^dx  >  /_  PT  -  x  I^Pl^dx  . 

-  ■^Fnn  4'n'  •'F  4'' 


Observe  now  that  P^^,  (tor  1  <  i  <  A)  are  uniformly  bounded  on  F  and  converge 

unlfomily  to  p,  (p’’.  Therefore,  we  just  have  to  prove 

lim  /p  Q({f,70^)dx  >  /p  Q((r,V0)dx  . 
n 

To  this  end,  we  write  ^  *  '*’n  obtain 

Q(0,V0  )  =  0(0,70)  +  0(0,71)1  >  +  20(0,70,71)1  ) 
n  n  n 

I  & 

where  Q(z,.,.)  is  the  symmetric  bilinear  form  associated  with  Q|z>  for  any  z  £  C  .  In 


particular  we  have 


0(0,70^)  >  0(0,70)  +  2Q(0, 70,71)1^) 


and  to  conclude  we  observe  that 

/p  O(0,'?0,’'l'„)ax  =  /  H*7i).^dx 

2  2 

where  H  is  some  fixed  function  in  L  while  7i))^  converges  weakly  in  L  to  0.  (A 
similar  proof  works  for  t,emma  III.3).  D 

We  next  consider  the  more  general  situation  of  a  nucleus  with  N  neutrons,  2 


protons  (so  A«N+Z).  We  may  number  the  wave  functions  in  such  a  way  that  9^, 


correspond  to  neutrons  while  ^ , . . .  correspond  to  protons.  We  also  denote  by  p^, 

T  ,  p  (x,y)  (resp.  p  ,  t  ,  p  (x,'/))  the  various  densities  of  neutrons  (resp.  protons) 
n  n  P  P  P 

1  •  G  • 

N  ”2  I’  « 

P  (x)  »  I  l».(x)l  ,  T  (x)  =  I  l7®  (x)|  ,  p  (x,y)  =  I  If  (x)9  (y) 

■’  1-1  ^  ”  i-1  ^  "  i-1  ^  ^ 


”2  *  *  2  ^ 

P  (x)  »  y  |9,(x)|  ,  T  (x)  -  I  |7(f  (x)l  ,  P  (x,y)  -  I  9.(x)if  (y) 
P  i=N+1  ^  P  i=N->-1  ^  P  i=N+1  ^  ^ 


In  this  general  casG^  the  functional  to  be  minimized  is 

2  t  X 

•28)  81^1 . Tj  =  /  ,  k  "  -  2“  ^  ■  ‘*0*  7>"’n  *  ^ 


•-V.Vi 
-•  V,*J 


9  i  OT  -  i  (P  T  >  0  T  )  +  ^  |7p|2  -  ^  (|7p  +  i7p„|^)  + 

4  4  n  n  p  p  16  16  n  P 

t  2  2 

9  ^  p  p^Ppdx  *  f-  /  /^  Op(x)  Pp(y)dxdy  -  |-  /  |  Pp(x,y)  l^dxdy 


*■  > 


a 

M 


-23- 


*•  '  e  R,  tg  >  0,  xg  e  <0,1),  a,  B,  y,  6,  t,  >  0. 

Observe  that  the  last  two  terns  obviously  correspond  to  a  Coulomblc  interaction 

between  protons.  Let  us  also  mention  that  vary  often  the  last  term  (Coulomb  exchange  term) 

is  neglacteds  this  makes  no  difference  on  the  type  of  mathematical  results  we  prove. 

Still  about  the  form  of  the  functional  it  la  worth  remarking  that  the  parameters  a,  B,  y,  B 

are  not  completely  Independent  since  in  fact  some  of  these  terms  are  exchange  tarns. 

The  first  situation  we  studied  in  this  section  corresponds  to  the  situation  when  the 

Coulomb  term  is  neglected  (e>0),  R“Z  and  P  “P  “TP»t  -  t  “tt.  Finally,  we 

n  p  A  n  p  ^ 

would  like  to  remark  that  , . . .  ,9^^]  is  not  invariant  anymore  under  all  unitary 

transforms  of  (e^,...,ej^)  but  only  under  the  transforms  of  the  form 

0  0 

«29)  0  -  („"  0  ) 

P 

where  (resp.  V^)  is  a  N  «  h  (reap.  Z  x  z)  unitary  matrix. 

And  we  consider  now  the  HP  minimization  problem 


(30) 


I  -  Inf{E[9^,...,?j^l  I  9^  e  h’(R-*),  /  j  pTdx  <  -  , 

R 

/  j  9j^»*dx  -  for  1  <  1,  j  <  N  and  for  M+1  <  1,  J  <  A} 


together  with  its  extension 


-  Inf{E(9,,...,9j^]  I  9^  «  h\r^),  /  ^  pTdx  <  -  , 

/  J  9j^9^dx  “  m^^  for  1  <  i,  j  <  H  and  for  N+1  <  1,  j  <  A) 


(31) 


H  0 

where  M  »  (m.  j)  *  (."  „  )  is  a  block  diagonal  hermitlan  nonnegative  matrix.  Using 

*  J  0  M 

p 

unitary  transforms  of  the  form  (29),  it  is  clear  that  we  may  still  diagonalize  M  into 


diag(m 


m. )  where  m,  >  0  for  all  i  e  (1 


a)  .  Therefore  I, 


Km 


m, )  where 


R 

^  3  “  ®l*ii  1  <  1,  j  <  N  and  for  N+1  <  i,  j  <  a} 

R  ^  ^ 

Observe  also  that  the  orthogonality  conditions  In  (30)  still  enable  us  to  write  down  Euler- 
Lagrange  equations  (the  HF  equations)  where,  up  to  a  unitary  transform  of  the  form  (29), 
the  matrix  of  Lagrange  multipliers  (e^  In  (17))  Is  diagonal  and  the  Lagrange  multipliers 
are  eigenvalues  of  self-adjoint  operators. 

Before  going  further  in  the  mathematical  analysis  of  (30),  we  would  like  to  mention 
the  way  the  parameters  tg,  t^,  Xg,  a,  6,  Y,  ^  are  chosen  In  realistic  computations.  The 
parameters  are  adjusted  by  a  simple  fit  to  the  binding  energies  and  equilibrium  densities 
of  some  fixed  nuclell  (essentially  oxygen-16  and  lead-25a).  Once  this  fitting  Is  performed 
(see  the  tables  in  D.  Vautherin  and  D.  H.  Brink  [48]),  one  can  compute  all  other  nuclei!  by 
solving  numerically  ( 30 ) . 

It  is  clear  that  conditions  on  the  parameters  are  needed  in  order  to  insure  that 
1  >  -•  (and  that  minimizing  sequences  are  bounded).  The  boundedness  of  I  and  the 
solution  of  (30)  are  analyzed  In  the 

Theorem  IV. 4;  1)  Assume  that  a  >  (B+6)/2,  <*  +  SAy  >  6+B.  Then  for  all  N,  Z  the  Inflmum 

I  e  (-<°,0]  and  for  all  R  >  0  there  exists  >  0  such  that  for  all  (<p^ , . . .  ,(^1^)  in 
the  minimizing  class  satisfying  E[e.| , . . .  ,9j^]  <  R  then  /  p+T  +  pT  dx  <  Cj^. 

11)  Assume  that  o  <  (B+5)/2  then,  for  all  N,  Z  >  1,  I  =  -®. 
ill)  Assume  that  the  conditions  given  In  1)  hold.  Then  any  minimizing  sequence  of  (30)  is 
relatively  compact, up  to  a  translation  if  and  only  If  the  following  condition  holds 
I  <  I(m^,...,m^)  +  I( 1-m^ , . . . , 1-m^>  for  all  e  [0,1)(1<i<A) 

(32)  ^ 

such  that  J  m  e  (0,A) 

1-1 

In  particular,  if  (32)  holds,  there  exists  a  minimum  of  (30). 


wywjrrjr^ 


Remarlcg ;  1)  The  analogues  of  the  remarks  given  after  Theorem  IV.  1  still  hold  here. 

11)  There  are  other  conditions  than  a  <  (B4d)/2  which  imply  that  I  -  We  mention 

only  this  one  to  emphasize  the  following  phenomenon!  take  M-Z,  e-0  then  In  this  case  It 


Is  often  assumed  In  the  Physics  literature  that  It  Is  enough  to  consider  (f  such 

that  P  “P  ”xP»''  without  changing  the  value  of  I.  And  this  Is 

n  p  2  n  p  2 

completely  false  In  general  (It  would  t>e  certainly  of  Interest  to  understand  completely 
this  kind  of  symmetry  breaking).  Indeed,  choose  a,  B,y,S>0  so  that 


a  <  (8+«)/2  ,  (a+Y)  >  (B+«)/2 


then  11)  Implies  that  I 


Inf{Ets^,...,ej^)  1  j  /  3  PT  < 


p  2 


T  <•  T  “  T-  T) 
n  p  2 


>  Inf{E'  1  9^  e  /  3  -  5^^,  /  3  pT  <  »} 


p2  +  I  |)px  +  (Y-  |)  |Vp  I  ^  + 


^-Ip^dx 

12 


Hence,  by  part  i)  of  Theorem  IV. 1,  the  restricted  Inflmum  Is  finite  as  soon  as 
(a+Y)  >  (8+5)/2. 

Proof  of  Theorem  IV. 4:  We  begin  by  proving  part  1).  We  first  observe  that  we  have 
denoting  by  Y'  =  yA6 

I  !  ‘Vn  ^  Vp’  "  ^  ^  " 

^  I  'Vp  *  Vn’  "  T-  'Vn  *  Vp’  *  ^  ‘”’n'’'’p’  '  '  I’^n ' '"1’^  l' >  ’ 

But  |Vp  1^<4PT  ,  |7p  |^<4PT,  and 
n  t\  n  p  P  P 

|(Vp^,VPp)|  <  2(P„t^  +  PpTp),  l(Vp„.VPp)|  <  2(PpT„  +p„Tp) 


■.  »  -  •U  •  e  •  •  **_fc.*^  *■ 


. 

.  - 


.■-'S'L'v'.'*" 


80  the  above  quantity  is  bounded  from  below  by 


( )(P„T  +  p  T  )  +  |•(P„T  +  p  T  )  -  ^  >  v  PT 

4  -  nn  pp  4np  pn  9  n  p 

where  v  is  a  positive  constant.  And  this  yields 

2  2 

]  +  ~  ^  ^  rJ'T  P  (*.y)^d*dy 


i  P 


for  some  C  >  0  (depending  only  on  tg,  Xq).  Next,  we  remark  that  pp(x,y)  <  Pp(x)pp(y) 
and  thus  by  standard  convolution  inequalities 


R^xR 


<  C(/  3  I  dx) 

R  i 


and  by  Sobolev  inequalitiea 


<  C(I  /  |Ve.l^dx)'^2  .  c(/  T  dx)'^2 


i  R 


where  C  denotes  various  constants  R^  depending  only  on  e  and  Ae  Therefore ,  one  gets 


EC9, 


,1  >  /  3  T  dx  t  V  /  PTdx  -  C  /  P^dx  -  C(/ 


T  dx/2 


It  just  remains  to  bound  conveniently  /  3  P^dx.  In  the  computations  that  follow  C 


denotes  various  constants  depending  only  on  A: 


/  j  p^dx  <  C  I  f  3|fj^|''dx 


i  R 


<  C(I  /  /  3ll’il’^'3x 


i  R' 


i  R 


by  Holder  inequalities,  and  since  /  -|t.|^<3x  =  1  for  all  i 


by  Sobolev  inequalities 


<  C(/  JVpI^ax)^'^ 

R 

<  C(/  3  PT  . 

R 

This  allows  us  to  conclude  the  proof  of  part  i) • 

The  proof  of  part  ill)  is  the  same  as  the  proof  of  Theorem  IV. 1  and  thus  we  will  skip 
it.  To  prove  part  il)  we  have  to  build  appropriate  test  functions.  By  simple 

considerations  it  is  enouqh  to  treat  the  case  so  denoting  by  f  ■■  i|i  ~  f2  ** 

2  2 

bave  ~  '  ^’p  ~  ^  since  we  will  take  real-valued  g,  i|i.  We  construct  spherically 

symmetric  functions  s”,  as  follows.  I.et  tg  >  1  be  such  that  -  «]tp  >  y- ,  let 

^0  ^  ^  ^0'  *^0  vill  have  to  be  determined  later  on  and  we  assume  at  least 

5q  <  1/4,  rg  <  1.  We  are  going  to  build  first  s*,  in  the  ball  B(0,rQ);  <i“  will  take 
values  in  the  interval  [tg-fig,tg]  and  f“  -  (2t^  -  compute  on  the  ball 

B(0,rg) 


s  '■  • !  »’'■  -  7  - 17 1  •  I'-pi'’  < 


where  F(t)  •  (2tg  -  ,  now  for  4g  small  enough  |f'1  takes  values  as  close  to  1 

as  we  wish  while  (e“)^,  take  values  arbitrarily  close  to  tg.  Therefore  fixing 

4g  >  0  small  enough  the  quantity  between  brackets  is  bounded  by  -v  with  v  >  0.  Since 
we  will  extend  y",  i(i“'  outside  B(0,tg)  in  such  a  way  that  e",  ♦“  and  their  first 
derivatives  are  bounded  by  fixed  constants  (depending  only  on  tg)  and  have  compact 
support  say  in  B(0,1)  we  deduce 

Els”,*®]  <  -V  • 

2  ^"’^0 

Now,  we  choose  tn  by  Imposing  {t„  +  F(t„-5.)  )  — —  <  1/8,  and  we  define  S®  as 


^  ^  <r*. *•*"•**•“  "-* '  •«*  *  ^  *••  *■“  *  ^j*  ’’  * "  “"’  1** 


follows  on  B(O,r0) 


♦■(X)  -  tg  -  «g  if  l*l<2^-S.  -t0  if 


t„  -  Sfl/S  +  j  m(|x|  -  if 

0  0  0  '  '  2  2  2n  '  '  2  2n 


r.  v.  />  V  •■ 'if ',V  "I'l*  VTWLTIVSIL^'^"t^VS!.V| 

Wi 

X* 

li-yw 


It  is  of  course  easy  to  extend  s"*,  outside  B(0,rQ)  as  we  claimed  edx>ve  and  we  can 


even  do  so  imposing 


/  j(»“)^dx  “  1,  /  (i(i“)^dx  «  1,  /  jS^^dx  -  0 

K  R"*  R'* 


(this  is  urtiere  we  use  the  restriction  on  tQ).  Coaq^uting 

..«|2 


2  2  4«  ,Xo  .  1  .3  ,'o 


1  .3, 


B(o,r„)  I’*  I  *«  -  V  r  "r  *  ill*  -  ‘r  -  i;^>  U  ^- 


we  prove  that  Bif”,  • 

ReB>arka»  i)  Let  us  observe  that  even  if  we  restrict  in  ii)  the  infimum  to  spherically 
symmetric  functions <  the  infimum  is 

ii)  The  idea  of  the  above  tedious  constructions  is  to  choose  at  least  locally  near  0 
P  ■  9*  +  constant,  jVfl  ~  |Vi(i|  ,  f  ~  4i  ,  thus  cancelling  the  iI^pI  term  while 
making  the  other  terms  a  pT  -  8(Pj,i^  +  ’’p^p*  ~  4  '  hpproxinately  equal  to 

(a  -  f  -  |■)PT. 

We  would  like  to  conclude  this  section  by  emphasising  Remark  ii)  following  Theorem 

IV. 4:  we  showed  there  that  if  N  •>  Z,  for  the  above  class  of  problems,  it  is  a  priori  not 

correct  in  general  to  restrict  the  infimum  to  configurations  such  that  t  =  t  “  t, 

n  p  2 

Pj^  “  Pp  ”  -j  P.  In  fact,  we  gave  an  example  of  a  dramatic  symmetry  breaking  in  the  Isospin 
variable  (between  n  and  p  i.e.  between  neutrons  and  protons).  Of  course,  this 
phenomenon  is  by  no  means  restricted  to  Skyrme's  Interactions  but  will  be  present  for 
general  interactions  (at  least  for  some  range  of  parameters).  \  precise  study  of  this 
phenomenon  (maybe  on  simple  model  problems)  certainly  remains  to  be  made,  investigating  in 
particular  the  possible  bifurcations  corresponding  to  it.  From  the  view  point  of  physics, 
this  symmetry  breaking  does  not  seem  to  have  been  observed  for  nuclei  such  that  N  =  Z  in 
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V.  other  lnt«ractlon> 

Ha  hava  consldarad  In  tha  pracadlng  aactlons  tha  casa  of  tha  Skyrma'a  Interaction 
which  basically  corresponds  to  tha  choice  (5)  (In  fact  tha  difference  between  neutrons  and 
protons  wave  functions  In  the  preceding  section  cosms  frca  the  fact  that  vave 

functions  should  depend  on  tha  so-called  Isoapln  variable  which  takes  fixed  different 
values  for  neutrons  or  protons  and  that  the  Hamiltonian  acts  also  on  this  Isospln 
variable).  It  Is  easy  to  understand  that  (5)  Is  a  very  simplistic  model  of  nuclear 
Interactions  which,  even  If  they  have  a  short  range,  do  not  have  "zero-range”.  However 
this  model  Is  often  used  because  It  already  gives  reasonably  good  numerical  results  and  the 
HF  equations  being  completely  local  are  somewhat  easier  to  confute.  Nevertheless  the 
theories  allowing  to  derive  the  effective  Interaction  V  from  first  principle  reveal  that 
one  has  to  consider  sure  sophisticated  parametrlzatlon  (l.e.  different  V...).  Furthepu>re 
several  extensions  of  the  HF  theory  (as  for  example  time  dependent  Hartree-Fock  problems, 
Rartree-Fock-Bogolyubov  theory. . . )  make  necessary  the  use  of  more  realistic  Interactions 
V.  On  all  these  l>aslc  Issues,  we  refer  to  J.  Kegele  [43]. 

In  this  section  we  are  mainly  Interested  In  the  case  when  V  Is  given  by  (3)  even  If 
It  Is  quite  clear  that  most  of  the  arguments  we  present  below  are  still  valid  for  much  more 
general  V  (including  (4)  as  another  example).  He  will  not  bother  to  Indicate  precisely 
what  are  the  mathematical  assumptions  we  need:  It  Is  an  easy  exercise  to  figure  out  In 
which  class  (for  Instance)  one  has  to  take  V  and  we  leave  It  to  the  reader. 

Again  to  explain  the  difficulties  encountered,  we  begin  with  the  scalar  case  which 
more  or  less  corresponds  to  the  case  of  the  alpha  particle.  For  V  given  by  (3),  we 
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/  jUrdx  -  X} 


whar*  X  is  a  positive  paramatex.  In  fact  solving  (35)  for  sons  X  >  0  asiounts  to  solve 
(34)  for  different  values  of  tg,  a,  S  as  it  is  easily  seen  by  a  scaling  arguaent*  We 

begin  by  a  simple  observation 

Proposition  V. 1 i  i)  For  all  R  <  «,  there  exists  <  »  such  that  Igl  .  ,  <  C-  if 

H  ^  ( R"* ) 

v  e  h’(*^),  B[e]  <  R,  /  ,Ul^dx  <  R. 

R 

ii)  One  has  always  <  0.  If  a,  6  >  0,  ■  0  for  all  X  >  0  and  there  is  no  niniaua 

of  (35). 

iii)  There  exists  Xg  e  (0,+«]  such  that  •  0  for  X  <  Xg,  Ij^  <  0  for  X  >  Xg. 

U  V 

Proof  I  Part  i)  is  easy  and  we  skip  it.  If  a,  6  >  0,  Bte)  >  0  for  all  s  0  and  thus 
aAnitting  I;^  <  0  the  remainder  of  ii)  is  clear.  To  prove  that  Ij^  (  0,  we  consider 
S  £  0(R^)  such  that  /  j|e|^dx  «  X  and  we  denote  by  “  e('^)o  Computing 


B[s-]  we  find 


\  ^  jlVel^dx}  ♦  i-  (tjj  ;  + 

C  R  C  R 

*7  f  f  Ul^<>‘)  V(a(x-y))|el^(y)dxdy 

^  <a  •* 


hence  Ete^l  *0  as  o  ♦ 

To  prove  iii),  one  first  remarks  that  1^^  is  nonincreasing  since  by  the 
concentration-compactness  argument  one  has  always 
(36)  I;^  <  1^  +  Ix_Y  for  all  y  €  (0,X) 

and  Ix_-y  *  0.  Therefore  one  just  has  to  prove  that  I^  “  0  for  X  small  enough.  Indeed 

if  e  £  /  ,|*|^dx  -  X 

R 

I  /  /  |p|^(x)V(x-y)  |9(^(y)dxdy|  <  Ifl^  lipl^  IVI 

L®(R^)  L^'^R^ 

<  CX  /  JVsI^dx 


r^jtrw 


'  •  j>  *.m. 


so  for  X  small  enough  E[e]  >  0  and  Hi)  is  proved. 

nie  proof  of  iv)  relies  on  the  following  choice  of  st  take  9  e  P(S^)  and  set 


9„(x)  -  a 
o 


-V4.,x, 


9(*).  'Rien  computing  we  find 


E[9g]  -  a 


-’/^(|;^/3|V9|^dx>.a-V^t„/3|9|^^^/^dx}. 

R  R 


*i  f  f  (o^V(o(x-y)))dxdy  . 


r^r3 


Remark  that  a^(ox)  ♦  (  '  +  ■-^■)ll  *q(x)  in  PMR^)  as  o  +  “  and  we  conclude 

p"'  |i"' 

easily  letting  a  go  to  +“.  • 

We  next  give  a  result  concerning  the  solution  of  (35) 

Theorem  V.2;  i)  Every  minimizing  sequence  of  (35)  is  relatively  compact  in  H^(R^)  up  to 
a  translation  if  and  only  if  the  following  condition  holds 

(37)  <  1^  +  I;^_Y  ,  for  all  y  e  (0,X)  . 

In  particular  if  (37)  )iolds  there  exists  a  minimum  of  (35). 

ii)  If  a,6  <  0  and  if  the  following  condition  holds 

(38)  Ij^  <  1^,  for  all  y  e  (0,X) 

then  there  exists  a  minimum  of  (35)  which  is  spherically  symmetric,  nonnegative,  smooth  and 
decreasing  with  respect  to  |x| . 

2 

Remarks!  i)  Very  little  is  known  on  the  values  of  X  (or  eauivalently  a,  6,  tn,-^*.>) 

2in 

for  which  (37)  or  (38)  holds.  We  only  got  very  partial  results  on  this  important  question, 
ii)  If  for  some  Xg  >  0,  (38)  holds  and  (37)  does  not  hold  then  there  exists  a  minimum  of 

(35)  while  some  minimizing  sequences  are  not  relatively  compact  even  up  to  a  translation. 

If  this  were  to  happen  this  would  be  an  extremely  interesting  situation. 

Proof  of  Theorem  V.2!  Part  i)  is  proved  by  a  direct  application  of  the  concentration- 
compactness  arguments  [33].  To  prove  part  ii)  we  first  observe  that  by  a  somewhat  standard 
symmetrization  argument  (as  in  E.  H.  Lieb  [28],  H.  Berestycki  and  P.  L.  Lions  [5])  one  sees 
that  I;^  agrees  with  the  Infimum  of  8(9)  for  9  e  H^(R^),  9  spherically  symmetric, 

nonincreasing  with  respect  to  |xl,  nonnegative  and  /  ^Ivl^dx  =  1.  Therefore,  it  is 

R 
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•nough  to  provo  that  if  (38)  holds  than  thare  axists  a  nlnimum  of 

if  -  lnf{E(f]  /  ?  e  f  la  spherically  synaatric,  /  |f  |^dx  “X},  the  other 

properties  of  the  miniaum  following  easily*  Mow  to  solve  if  we  may  either  apply  the 

concentratlon-coe^actness  arguments  In  presence  of  symetrles  (see  [34] )  and  conclude 

observing  that  since  vanishes  for  X  small  then  11m  n  “  0  for  dlX  X  >  0  and 

n 

thus  (38)  Is  equivalent  to  if  -  <  l"  +  11m  n  ^(X-Y)/n*  **  *  “ore  standard 

n 

line  of  arguments  showing  first  (as  In  [28]*  [46]*  [S] )  that  there  exists  a  minimum  Sg  of 

Inf{E[9]  /  <p  €  H^(R^)*  9  Is  spherically  syonetrlc*  /  lej^dx  <  X)  and  concluding  that 

R^ 

/  ^  Since  (38)  holds* 

R 

Before  going  into  the  general  case*  we  study  problems  lllce  (3S)  with  ft  ••  0  In  which 
case  (35)  reduces  to 

(39)  I.  -  Inf{t  /  |9|^‘^^''^dx  +  4  /  /  |•|^(x)V(x-y)|9|*(y)dxdy 

“  R**  ^33 

9  «  L*(H?)  ^  l'^'^’(B?)*  /  kl^dx  -  X> 

14 

(in  fact  the  value  j  plays  no  role  in  the  analysis  below,  for  exanple  everything  below 

14 

remains  true  If  we  replace  —  by  any  p  >  4)*  We  still  denote  by  E[9]  the  functional 

that  we  wish  to  minimize.  We  can  prove  the 

Theorem  V.3!  1)  If  for  any  P€lVr^),P>0  vm  have 

(40)  f  /  p (x)V(x-y)p (y)dxdy  >  0 


then  1;^  >  0  for  all  X  >  0  and  there  Is  no  minimum  of  (39)*  On  the  other  hand*  If  (40) 

does  not  hold  for  some  p  c  L^(R^)*  P  >  0  then  Ix  <  0  for  all  X  >  0*  This  Is  the  case 

(X  B 

if,  for  example,  a+B  <  0  or  - 7-  +  — rr  <  0*  In  all  that  follows  wo  assume  that  Ix  <  0 

N/2  N/2  * 

U  P 

for  all  X  >  0. 

11)  Every  minimizing  sequence  of  (39)  is  bounded  in  L^(R^)  O  L^^^^(R^). 

ill)  Every  minimizing  sequence  of  (39)  Is  relatively  conqpact  in  L^(R^)  H  l'^^^(R^)  up  to 

a  translation  if  and  only  if  (37)  holds.  In  particular  if  (37)  holds  there  Is  a  minimum  of 


iv)  The  condition  (37)  holds  if  X  is  small  enough 


v)  If  a  and  8  are  negative ,  (37)  holds  for  all  X  >  0  and  there  la  a  minimum  of  (39) 
which  la  spherically  symmetric,  nonnegative,  nonlncreaalng  with  respect  to  |x|  and  with 
compact  support. 

Remarlcs  t  1)  If  wo  assume  that  (37)  holda  at  X  >  0  than  we  can  prove  that  either  (37) 
holds  In  a  neighborhood  of  X  or  there  exists  a  minimum  of  (39)  -  that  wo  denote  by  s-  - 


that  »o  ^  ^  l''^/^(*^)  and 

14  11 

E’(»o>  “  0  !■•••  ~  ^0  ’O 


a. a.  In 


(we  may  always  assume  that  Vq  la  real-valued,  nonnegative). 

Indeed,  If  for  any  sdnlmum  of  (39),  (41)  does  not  hold  than,  assuming  that  we  have 
built  a  seguence  J  X  such  that  ( 37 )  holds  for  1^  and  denoting  by  the 

associated  minima,  on  one  hand  converges  (up  to  subseguences  In  j^14/3(,3)  ^ 

to  some  minimum  ipg  of  I^  and  on  the  other  hand  there  exist  6^  satisfying 
(42)  *  ^n’n  "  °  R^,  0  <  v  <  6^  <  C 

for  some  positive  constants  v,  C. 

Row  we  argue  by  contradiction:  if  (37)  does  not  hold  for  where  then 

n  n 

there  exist  7^  >  Xj^/2  such  that 

^X  ■  S  ^  ^X  -Y  '  °  <  ^n  <  ^ 

n  n  n  n 

and  alnca  (37)  holds  for  we  have  that  ♦  Xe  Hext^  If  (37)  does  not  hold  for  1,^  , 

there  would  exist  €  (Yn/2'^n'  *  ^y  -6  *  particular  we  have 

n  'n  n 

Ia  ■  ^  *  "x  -Y  ■  *  S  ^  "x  -Y  • 

n  n  n  n  n  nn  nn 


But  we  always  have 


_s  ^x  -r  *  ^x  -5 
n  n  n  n  n  n 


so  the  above  eguallty  yields 


*  h  -«  '  ^X  -6  “  s  -«  ^  ^X  -Y  • 

n  nn  nn  nn  nn 


Since  (37)  holds  for  Ii,  the  first  equality  Ijnplies  that  5  ♦  X-  But  then  the  second 

A  n  ^ 

equality  gives  a  contradiction  since  (37)  holds  for  X  small  enoughs  Therefore  (37)  holds 


To  conclude  we  argue  as  in  [34] :  observe  that 


n  n  n 


*“  -‘  n  “ 


<  I  -  K(X  -Y„) 
Y  n  n 

n 


for  some  K  >  0  (the  almost  equal  sign  can  be  easily  justified  and  the  above  Inequality 
holds  rigorously).  On  the  other  hand,  we  prove  below  that  converges  to  a  negative 

constant  as  X  goes  to  0.  Hence  we  get 


<  C(X  -Y„) 
n  n  n  n 


and  the  contradiction  proves  our  claim.  • 

2)  We  would  ll]ce  to  remark  also  that.  If  we  do  not  assume  that  a  and  B  are  negative, 
the  question  of  the  spherical  symnetry  of  the  minimum  of  (39)  (when  It  exists,  as  for 
example  when  X  Is  small)  Is  open. 

3)  Let  us  finally  point  out  that  somewhat  related  problems  are  considered  In  J.  F.  G. 
Auchmuty  and  R.  Beals  [1],  [2],  P.  L.  Lions  [36],  [33]. 

Proof  of  Titeorem  V.3!  The  proofs  of  1)  and  11)  are  standard:  the  sign  question  being  a 

consequence  of  the  difference  of  homogeneity  of  the  two  terms  In  E,  and  the  negativity 
•N/2  •N/2 

of  when  au  +  Bv  <0  being  proved  as  In  Proposition  V.  1.  Part  ill)  Is  proved 

by  a  simple  application  of  the  concentration-compactness  method. 

We  now  prove  part  Iv).  We  first  show  that 

(43)  I^x"^  I  -  Inf{l  /  /  |e|^(x)V(x-y)|9|2(y)dxdy  /  9  e  L^(R^),  /  jUl^dx  -  1}  . 

* 

Observe  by  the  way  that  the  Inflmum  In  I  Is  achieved  by  a  simple  application  of  the 
concentration-compactness  principle.  To  prove  (43)  It  Is  enough  to  remark  that  on  one  hand 
Ij^  >  X^  I  while  on  the  other  choosing  9^^  c  p(R^)  such  that 

^  ”  T  /  /  l'P_l^(x)V(x-y)|9  |^(y)dxdy  <  I  +  ^ 

^  n  ^  li  r#  n 

R^xr^ 


and  f  3 1  I  ^dx  =  1 , 


we  obtain 


I.  <  E[/X  9  1  <  C  +  X^(I  +  -) 


Next,  let  ua  argue  by  contradiction  to  prove  part  iv):  assume  there  exist  + 


and  Y„  e  such  that 


This  yields 


n  n  n  n 


^X  ^Y  Y  2  *x  -Y  X  -Y  2 

—S.  .  -In.  r-iX  +  "  "  f  »  "1 

2  2  '■X  ^  X  -Y  '  X  ' 

X  Y  n  n  n  n 

"  "  ^n 

and  this  combined  with  (43)  istpliea  that  up  to  subsequences  y"  converges  either  to  1  or 

n 

to  0.  Replacing  if  necessary  Y^  by  Xj^-y„  we  may  always  assume  that  Y^/X^  converges 
to  1.  Next,  we  fix  n  and  we  talte  a  minimising  sequence  for  1^^.  We  remark  that 


A  /5  A  ^ 

<  ii«*[(:^)  e^l  *,1%!  +  (— )  =2^^^ 


n  k  n 


where  we  denote  by  S^tsl  “  tg  /  jlfl  ®h  the  other  hand  [9j^)  +  E2^^k^  ^  ^Y^ 


so  we  obtain 


X  2+1/3  X^  2  X  2+1/3 

Ix  <  S  -  u»  l*2tv(^r)  - 

n  'n  'n  k  n  n 


X  2  X  1/3 
ni  rr  n-, 


<  1,^  +  ((-^l  -  ■ 


Recalling  (44)  we  deduce  finally 


X  2  X  -Y  1/3 
r  n>  If,  .  n  n^v 


..  <  ((tt^  -  +  (r)  Hi  +  -“)  -  ■>}  iiB  -®2'V^  • 


Of  course  lim  {l  -E_(e.  )}  depends  on  Y_  and  similar  arguments  to  those  used  to  prove 
k  ^n  ^ 


(43)  show  that 


iia  -Ej(e^)}  y;  +  0  . 

R  n 

^X  -Y  ^x  -Y  X  -Y 

Next  since  77-^^ — 7 — -  - ^  ^  °  dividing  the  above  inequality  by  (X^-Y^lY^  «' 

'V^n’^n  (X„-Y„)  "'ll  " 

n  n 

obtain  passing  to  the  limit i  0  <  21,  contradicting  the  negativity  of  I.  And  the 


contradiction  proves  iv) . 


We  now  conclude  the  proof  of  nieorem  V.3  by  proving  part  v).  The  proof  Involves 
several  steps :  1)  we  show  the  existence  of  a  spherically  sysnetrlc ,  nonincreasing  minimum 

of  (37),  2)  that  such  minima  have  compact  support.  It  Is  easy  to  conclude  that  (37)  holds 
by  observing  that  if  Vj  are  the  respective  minima  of  1^,  for  some  y  e  (0,1) 

and  If  say  93  supported  In  a  ball  of  radius  R  then  considering 

9(x)  -  e^(x)  +  ^^(x  +  2Re) 
where  e  Is  any  unit  vector,  we  obtain 

/  ^ 

R  R  R 

E(el  «  E[(P^1  +  +2  ^  ^  (x)V(x-y)e2(y+2Re)dxdy 


and  (37)  Is  proved.  (Observe  that  ^  0  >nd  V  <  0.) 

To  prove  1 )  we  argue  as  In  the  proof  of  Theorem  V. 2  Introducing  the  problem  In 

p  “  lei ^ 

=  Inf{tg  /  jlpl’^^^dx  *  ^  f  !  P(x)V(x-y)p(y)dxdy  / 


R 

.1/-3, 


r^r3 


p  spherically  symmetric,  p  €  l'(R^)  O  p  >  0  a.e.,  /  ,pdx  <  X}.  Then  this 

vr 

problem  Is  solved  exactly  as  In  P.  L.  Lions  [36]  using  the  spherical  symmetry  and  the 


smoothing  properties  of  the  kernel  V(x-y)  and  there  exists  a  minimum  Pq  which  is 


nonincreasing  (using  again  symmetrlzatlon  arguments).  If  we  prove  that  /  ^  Pq^  **  ^  then 


Step  1 )  Is  completed  considering  p 


''Pjj.  In  order  to  do  so  we  argue  by  contradiction 


and  we  assume  that  /  ^  Pg^  ^  Then  the  necessary  conditions  for  minimality  may  be 


written  as 


y  tg  Pg^^  +  j(pg*V)  =  0  s.s.  oo  ths  set  (Pg  >  0} 


y  tg  Pg'^^  +  ^  °  ^**0  “ 


But  Pg  Is  spherically  symmetric,  nonincreasing  with  respect  to 


so  the  set  {pn  >  0} 
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'"iw  ■  ^ 


-  |Pp(x,y) |^}dxdy  +  ^  p^^^dx  . 


whera  W^,  H^,  Nj^,  (1-1,2)  are  given  constanfca  (which  in  practice  are  not  independent  - 
roughly  speaking  »  -B^,  -  -M^)  and  V^(i-1,2)  are  given  byt  Vj^(x)  -  exp(-|x| ^/ll|) , 

and  u^,  U2  ^  constants.  Of  course,  we  are  using  the  same  notations 

concerning  p,  p^,  p^,  t  as  in  the  functional  (28)  for  Skyrme's  interaction  (section  IV). 
nie  HF  minimization  problem  is  then 

(46)  I  -  Inf{Et9^,...,Sj^l  /  <  H’(m^),  /  j  ^i’j***  “  ®ij 

for  1  <  i,  j  <  N  and  for  N+l  <  i,  J  <  A}. 

We  %fill  not  state  a  result  on  this  problem  because  exactly  the  same  result  as  in  part 
iii)  of  Theorem  IV. 4  holds  here  (and  the  remarks  following  Theorem  IV. 1  or  Theorem  IV. 4 
also  hold  here).  Of  course  since  V^  e  l‘‘(K^)  for  i  -  1,2  t)ie  infimum  and  minimizing 
sequences  are  automatically  bounded  and  contrarily  to  Theorem  IV. 4  no  restrictions  on  the 
coefficients  need  to  be  made  prior  to  the  analysis  of  minimizing  sequences. 


> 

-•  j" 


1 
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VI.  Spin-orbit  forces 


XJp  to  now  we  heve  constantly  ignored  the  spin  dependence  of  the  various  wave 
functions.  However,  if  this  omission  greatly  simplifies  the  presentation  and  (probably) 
the  mathematicF  f  the  HF  minimization  problems,  for  practical  and  realistic  computations 
one  has  to  cope  with  the  spin  dependence  and  its  consequences:  the  spin-orbit  force.  It 
is  our  goal  here  to  try  to  explain  the  form  of  the  spin-orbit  force  and  to  show  that  in 
order  to  have  a  bounded  infimum  some  precautions  have  to  be  talcen  (and  this  does  not  seem 
to  have  been  always  the  rase  in  the  Physics  literature  on  this  matter).  Let  us  also 
mention  that  in  the  remaining  sections  of  this  paper,  we  will  again  slip  the  spin 

dependence  even  if  it  can  loe  restored  without  affecting  the  mathematical  results  (provided 

t 

one  considers  spin-orbit  forces  with  the  appropriate  restrictions  described  below) . 

First,  we  explain  how  wave  functions  depend  on  spin  and  we  will  do  so  by  only 
explaining  the  computational  rules.  In  everything  we  said  in  section  II,  one  has  to 
understand  now  that  *(x^, . .  .,Xj^),  9^(x)  in  fact  depend  on  other  variables  namely 

♦(x^,o^»  X2,02;...»x^,Cj^) ,  Sj^(x,(j) 

where  the  spin  variables  ta)ce  only  two  values  say  +1  and  -1.  If  we  denote  by  x^  * 

<  i  <  A),  the  Pauli  principle  now  states  that  the  antisymmetry  condition  (2)  has 
to  loe  understood  now  as  a  condition  on  permutations  of  the  variables  x^(1<i«A).  Then  the 
remainder  of  the  derivation  of  HF  problems  goes  through  as  before.  It  is  possible  to 
consider  now  s^(x,o)  as  a  pair  of  complex-valued  functions  (spinor)  that  we  will 
indifferently  denote  by  (»^(1),  ?^(-1))  or  (et,9T). 

The  orthogonality  condition  becomes  (if  no  differences  between  neutrons  and  protons 


are  made) 


/  j  e^(1)9j('')  +  9^(-1  )9j  (-l)dx  “  6^^,  for  1  <  1,  j  < 


The  spin  dependence  affects  the  Hamiltonian  H  and  the  potential  V  in  two  ways:  the 
first  one  is  through  the  so-called  spin-exchange  operators  (P^)  which  will  basically  mix 
the  various  products  of  The  second  one  is  more  dramatic;  it  is  the  so-called  spin- 

orbit  force  which  can  be  thought  of  as  an  additional  two-body  term. 
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Typical  aodala  of  tha  spin-orbit  fores  ars  ssro-rangs  Modsls  ooniparabXs  to  ths  choice 


(5)  of  potentials*  This  sodel  leads  to  a  functional  which  in  the  case  of 

Shyrme's  interaction  (see  section  IV)  is  given  by 


Stv. 


-*'’1 . V 


(48) 


+  ^  /  ♦  (Vp^.J  )  +  (Vp^.J  )dx 

2  ^3  n  n  P  P 


where  E  1  given  by  (28),  the  parameter  8  already  occurs  in  B,  Nq  Is  a  positive 
pareuneter,  J  is  the  so-called  spin  density  that  we  describe  below  and  3^,  3^  are  the 
spin-densities  for  the  neutrons  and  protons  and  are  built  in  the  same  way  as  J 
restricting  the  various  sums  to  the  neutrons  or  protons  wave  functions  9^  (as  we  did 


for  p^,  Pp. ..).  Let  us  also  mention  that  the  densities  T,  p  now  mean  of  course 


1-1  ^  1-1 


We  now  describe  Js  J  is  a  function  talcing  values  in  iP  which  may  be  written 

as  (see  [48]) 

(49)  J(x)  -  (-i)  I  9!(x,o)IV9^{x,o')  *  <o|o|c'>l 

j,o,o*  J  ^ 


where  a  is  the  Pauli  spin  matrix.  The  above  bracket  means  that  <o|$|o'>  is  a  point 


of  whose  coordinates  are  the  results  of  the  action  of  the  2x2  matrices  a^,  a^, 


described  below  on  (a, o')  where  the  spin  variables  o,  o'  take  values  now  in  (.)  or 


(°)  with  the  conventions  "+  ”  (g*"'  ""  ”  ^i*"’  matrices  o^,  a^,  are  given  by 


(50) 


0  *  ^  “i  0*  *  ^0  *1  " 


For  example  if  o  -  o'  -  (  )  then 


<o|a^|o'>  -  (1  0)(®  ’)(’)  -  0 

<o|o^|o'>  -  (1  0)(_J  -  0 

cojo Jo  -  (1  0)(’  °){’)  -  1 


In  general  the  above  quantity  J  ia  not  real  and  at  thia  point  we  need  to  explain  an 
important  assumption  in  Hartree-?ock  theory  in  Nuclear  niysics*  Mhen  including  spin 
dependence  and  spin-orbit  forces  (as  one  should),  one  has  to  work  in  the  case  of  even-even 
nucleii  i.e.  one  assumes  that  N,  Z  and  so  A  are  even.  Futhermore  one  assustas  that  the 
subspace  of  occupied  single-particle  states  is  invariant  under  time  reversal  and  this  means 
mathematically  that  for  all  j  (j  «  or  j  «  {N+1,...,A>)  there  exists  j'  (J*  e 

or  j*  «  {N+1,...,a})  such  that 

(51)  Sji(x,o)  -  -o  Sj(x,-a)  for  all  x  e  R^,  o  -  jfl 

+  ■*  «  +* 
i.e.  s j I  ■  "Fj  *  e j I  “  S j  • 

It  ia  possible  to  use  thia  assumption  by  dividing  by  two  the  number  of  unknowns 

(N,  Z,  A  become  N/2,  Z/2,  A/2)  and  we  still  denote  by  N,  Z,  A  those  reduced  numbers)! 

then  the  HF  minimisation  problem  and  the  functional  remain  the  same  and  one  may  compute  the 

three  components  J„,  J.  of  J.  A  tedious  computation  yields 

X  y  s 


-  X  {lm(«|.^(1)7  »  (D)  -  lm(F^(-1)V  e^(-1))  + 


+  Re(e^(1)V^9^(-1))  -  Re(s^(-1)V^S^(1))} 


J  -  I  {-lm(9*(1)V  s,(1))  ♦  lm(s,(-1)V  .,(-1))  + 
y,^.  1  xi  1  <1 


+  Im(s*(1)V  v,(-1))  +  Im(s,(-1>7  ..(I))} 

1  Z  1  L  Z  X 


(54) 


*  *  • 

J  -  I  {-R*(»,(1)7  *,(-1))  +  R*(f,(-1)V  f(1)) 

X  XX 


«•  will  only  Invoatlgate  hera  the  dlCflcultlas  concerning  the  boundednees  of 
mlnlmlElng  sequences  and  the  flnlteness  of  the  InflBum,  casing  frost  the  addition  of  the  two 
spin-orbit  terms  In  the  functional  E  given  by  (48).  Of  course,  we  are  Interested  In 


1  «  inf{Ete^,...,ej^]  1  *^  €  h\*^),  /  jpxdx  <  «, 

R 

I  /  3  e^(x,o)e^(x,o)dx  -  for  Kl,  J<tl  and  for  M+1<1,  J<h} 


Once  these  equations  are  solved  positively,  then  the  analysis  of  (55)  goes  along  the  sasM 
lines  than  In  the  preceding  sections  (and  raises  even  more  open  problems). 

The  considerations  we  give  below  show  that  I  Is  finite  In  the  case  of  the  S)cymw’s 
Interaction  If  Wq  Is  small  enough  (and  a,  S,  5,  Y  satisfy  the  conditions  of  1)  In 
Theorem  IV. 4),  while  for  other  Interactions  having  finite  ranges  (but  no  sero)  I  Is  never 
finite.  The  conclusion  Is  that  the  spin-orbit  force  cannot  be  ta]cen  as  a  sero-range  two- 
body  Interaction  and  one  has  to  use  Instead  spin-orbit  force  term  ll)ce 

-(r,-r3)V  ,  , 

e  [(r^-r^)  *  (v^  -  v^)  •  (o^  +  o^)).  We  will  not  try  to  explain  to  non-expert 

readers  what  this  terms  means i  let  us  just  mention  that  it  leads  to  HF  minimisation 
problems  Involving  terms  ll)ce  the  ones  we  are  analysing  except  t)iat  these  terms  are 
nonlocal  and  so  present  no  more  singularities  nor  unbounded  features. 

By  Inspecting  the  proof  of  parts  1),  11)  of  Theorem  IV. 4,  one  chechs  easily  that  part 
11)  still  holds  for  E  and  thus  we  will  assume  that  a,  6,  6,  y  satisfy  the  conditions  of 
part  1)  of  Theorem  IV. 4.  Therefore  we  find  that  if  (93,...,*^^)  are  In  the  minimising  set 


".*.*'*’  -m*  .r-. 


>v/jT  +  PTdX-C+|  /  +  |jj 


♦  ^  ><** 
2  —3  n  n  P  P 


for  soao  constants  v,  C  >  0*  If  follows  sasily  that  if  Wq  is  soall  (Wq  <  8vB),  I  is 

finite  and  If  <  R  than  /  ^  t  -f  pt  dx  <  for  some  poattlve  constant 

R 

How,  we  are  going  to  show  by  an  exaag>le  that  I  la  no  more  finite  (l.e.  I  ~  -*‘)  If 
we  consider  more  realistic  Interactions  such  as  the  ones  considered  In  section  V.  To  be 
more  specific,  we  consider  the  functional 

W 

(56)  Bte,,...,ej^]  -  ♦  2  / 

where  B  la  given  by  (45).  He  claim  that  this  functional  is  not  bounded  from  below  on  the 
minlstlzlng  set  {f.  «  H^(H?),  I  /  ,  “  *i-i'  1  <  1,  j  <  H,  H+1  <  i»  j  <  M  •  In  fact, 

we  believe  that  related  exaa^les  show  that  even  In  the  case  when  E  is  given  by  (48) 

I  X  -<•  If  Wg  la  not  SBMll  enough.  To  prove  our  claim,  we  begin  by  a  simple  scaling 
argusMnt  where  X  >  0  Is  the  scaling  parameter,  (e^,...,f]^)  Is  any  test  function  in  the 
minimising  set 

. /  a  ’  «x]  *  c,  t 

X  R 

2  t 

*  T  r  { ■  I  ^  7?  * 

R^xR^ 

1  "o 

♦  -1^  /  j(Vp,J)  ♦  +  f’Pp'^p^*** 

X  R 

where  Is  a  bounded  constant  (depending  on  f^,...,9|^).  The  example  below  shows  that 

the  last  term  may  be  negative  and  thus  our  claim  Is  proved  sending  X  to  0.  As  In  the 
proof  of  Theorem  IV. 4  it  Is  enough  to  build  e^,...,e,^  near  0  and  we  will  actually  build 
e^,...,S^  so  that  T,  p,  p**  (for  q  <  9),  |Vp|*  are  integrable  near  0  while  p(dlv  J) 
has  constant  sign  and  Is  not  Integrable  at  0.  Then  it  Is  easy  to  approximate  and  obtain 
values  which  go  to  Our  choice  of  9^,...,V]^  is  the  following;  for  1  >  3  take 


."t  .-v  v-'.--  -• 


>  J*  /  #'kl 


to  be  0  near  0/  for  1  <  2  taXe  to  be  0  near  0  and  real  near  0*  Denote 

by  V  ••  92'  ~  V'  'I'v  *^2  r«al.  He  find  that 


W  W 

!  j(Vp,J)dx  “  “  2^  /  3  P  'U.v  J  dx 

R  R 


and  nearby  0,  dlv  J  reduces  to  2  i|(*  i|i).  Observe  also  that  the  term  in 

p  div  J  given  by 


f  e-i2  .  j2.^  c  .  3.  3.  A 

^  *2  3j;  ^  -  sr;  ♦2  air;  ♦i**'  “  ® 


2.  3 


Hence,  tre  only  have  to  look  at  T  -  9  (-jj—  <1^  near  0.  We  next 


choose  9(x)  -  jxj^*  near  0  for  some  a  <  0  to  be  determined  later  on,  and  we  taka 
t)i,(x)  -  ?(x)  cos  6(x),  ip2(x)  ”  5(x)  sin  e(x).  So  T(x)  -  |x^l“c(x) 

[j^  C  3^  0  -  0  3^  Cl*  W®  finally  choose  c(x)  -  (x^  +  bxj  +  Xj)®'^^,  9(x)  -  |x|^'^^ 


1  ""2  ’"1  2 

with  b  >  0,  by  1  and  the  exponents  6,  Y  will  be  determined  later  on.  With  these 


choices  f  Is  given  by  f  -  8Yr“8®x^X2( 1-b)s®“^r^”^  where  s  -  (x^  +  bx|  +  x|)'^^. 
Therefore  If  is  not  integrable  at  0  if  2a  +  26  +  Y  <  “1j  while  t,  p,  p**  (for  q  >  1) 


vp 


are  integrable  if  o  >  -1/4,  6  >  -1/4,  a  >  -3/(2q1,  6  >  -3/(2q),  B+Y  >  -1/2.  Then 


if  q  <  9,  choose  o,  6  <  -  ^  near  "  IT  <  "  "J  "““’f  ”  "J  then  all  the  above 


conditions  are  satisfied. 

We  would  like  to  conclude  this  section  by  Inspecting  the  size  of  the  spin  orbit  term 


—  /  Vp  .  J  dx  in  the  case  of  a  spherically  symmetric  configuration  (the  precise  meaning 


of  that  choice  will  be  given  in  section  VIII).  Following  Vautherin  and  Brink  [48]  wo  see 
that 


J  -  <2V”»a‘V’»  -  w’’ 

•  '  An  IT  « 


4xr  a 


while  p(x)  »  - -  \  (2j^t1 )R^( I x| ) ,  where  to  some  positive  integer,  j*  *  *a  t  V2 

4'rr  a 

and  the  sum  over  a  means  the  so-called  sum  over  occupied  states  (the  set  of  levels 


compatible  with  the  numbers  of  nucleons  -  see  section  VIII).  Hence,  the  spin  orbit  term 
gives 


2 


/  j  Vp.J  dx 


<x 


1 


•  (I  • 

6 


If  w*  chooa*  only  on*  occuplad  atat*  a  (thla  My  alwaya  ba  aehlavad  by  taking  varloua 


Rj,  with  dlatlnct  aupporta  in  (0,«))  with  la  “  *o  ■*■  Y  '  '** 


daduca 


^  /_3  VP.,  dx  .  ^  rXK  -  7  -T 


and  tha  abova  acallng  argument  ahowa  that  awan  In  tha  context  of  apharlcally  ayamatric 
conftguratlona,  tha  apln-orblt  term  la  "too  unbounded*.  In  fact,  tha  abova  computation 
alao  ahowa  that  even  for  Skyrma'a  Interaction  condltlona  on  Wq  (compared  to  tha  other 
conatanta  6,0...)  have  to  ba  Impoaad  In  order  to  have  a  meaningful  HF  mlnlmlaatton 


problem^ 


VII.  A-body  problema  in  tluclear  Phyaica  and  Biomas-Fermi  approximations 


In  this  section,  we  first  make  a  few  comments  on  the  translation-invariant  A-body 
problesM  of  the  form  (1)  and  on  the  minimization  problem  (S).  Then  we  investigate  the  role 
of  the  density-dependent  term  in  the  Hamiltonians  H  which  are  being  used  in  practice  in 
nuclear  Physics.  Finally,  we  conclude  this  section  by  discussion  the  validity  of  HF 
approxisMtion  and  we  briefly  discuss  the  Thomas-Fermi  approximation. 

We  begin  with  problem  (8)  where  H  is  given  by  (1).  In  the  reMrks  which  follow  we 
will  not  bother  given  precise  assumptions  on  V  which  guarantee  an  eas.  justification  of 
the  arguments  below  (again  it  is  an  easy  exercise  that  we  leave  to  the  reader).  We  first 
observe  that  (8)  has  never  a  minimum:  indeed,  let  t  e  H  with  /  |*|^dx  >1  we 

consider  for  X  >  0 

(57)  ♦(x,,...,x.)«  —  ♦(x,  -  x  +  -fx,...,x  -  x+  -fx) 

A  1  A  I  A  A  A 


where  *  ”  "T  I 

”  •  a 


*i’ 


One  checks  easily  that  e  ((  and  /  dx  -  1.  Hext  tra 


compute  E(*\)  *  2  obtain 


A  R  1 


Hence,  if  t  is  a  minimum  of  (8)  the  above  equality  implies  I  -  0  and  this  is  not 
possible  since  t  €  4/0. 

The  above  equality  also  shows  that  the  ground  state  energy  E  is  also  given  by 


(58) 


2  2 
rft  V  f  In  -1) 


"  =  ^  ^  ^.3aI?  Vl  ""  " 

i  R  R  X 


^  ^  eH  nH’(R^*),  /  jj^|*|^dx  «  1)  . 

i<j  R  ^  R 


And  now  the  translation  invariance  does  not  imply  anymore  a  priori  that  minima  do  not 
exist.  In  fact  to  our  knowledge  no  existence  results  of  minima  for  A  >  3  are  known  for 
the  above  problem.  Observe  also  that  the  above  quadratic  functional  is  invariant  under  the 
transformation  (4  4.)  for  X  >  0. 


In  fact,  tn  practice  it  auiy  ba  important  to  apply  HF  method  to  the  above  functional 
Instead  of  and  when  we  Inject  Slater  determinants  Into  the  above  functional  we 

obtain  the  following  quantity  (which  clearly  replaces  the  term  ^  /  Tdx  that  we  had  in 
the  preceding  sections) 

(59)  -i’ 

and  everything  we  did  In  the  preceding  sections  Is  easily  adapted  to  this  new  situation. 
We  now  uhe  two  ramarlcs  for  Improving  the  Physics  applications  of  HF  methods.  First 


of  all,  there  Is  a  slight  difference  of  mass  between  neutrons  and  protons  and  this  could  be 

* 

Incorporated  In  H  and  In  everything  we  did  before  by  replacing  -  I  &  by 
2  N  2  1-1  1 

-  q—  J  A  -  J  A  where  m_,  m  denote  respectively  the  masses  of  neutrons 

2IB  .  -  Xj  2fll  “  P 

n  1-1  1  p  i-M+1  1 

and  protons. 


\  A  where  m^,  m  denote  respectively  the  masses  of  neutrons 


The  next  resMr)c  concerns  the  density  dependent  terms:  in  the  preceding  sections  (III 
to  VI)  all  HF  functionals  Incorporated  terms  nonlinear  with  P  homogeneous  of  a  degree 
different  from  4  and  obviously  such  terms  cannot  be  obtained  from  through  the  method 
presented  In  Section  II.  In  fact,  to  Improve  the  numerical  computations  obtained  through 
HF  methods  Huclear  physicists  have  added  to  the  Hasdltonlan  H  phenomenological  terms  of  the 


(60)  t,  I  V(x.-x,)V(x,-Xy) 

l<j<)t  ^  ^ 

where  V  for  example  may  be  V  -  6-.  Now  if  we  use  Slater  determinants  this  terra  gives 


some  term  ll)(e 


r  b,  /  /  /  P(x)V(x-y)p(y)V{y-z)p(i)dxdydr 
V  3  ,  ,  , 


or  T  b-  /  -  P^dx  If  V  -  5g.  Recalling  that  we  are  suppressing  the  spin  dependence  one 
sees  that  such  a  term  Is  equivalent  on  Slater  determinants  to  a  two-body  density-dependent 


Interaction 


J  bj  p(  6(x^-x^) 


(all  this  Is  formal  because  the  absence  of  spin  does  rnahe  matters  a  bit  trivial).  Roughly 


spealclng  the  term  (60)  provides  a  simple  phenomenological  representation  of  many-body 


•ffacta  and  la  auppoaad  to  daacrlbe  the  Influence  of  all  other  nucleona  to  the  Interaction 


between  two  of  them, 
conalder  (aee  H 


It  haa  alao  been  obaerved  that  Inataad  of 
.  Bethe  [8] )  which  leada  to  the  term  / 


p  It  la  often  better  to 
p  2+2/3^  uaed  In  aectlon 


V.  Of  courae  a  term  like  (61)  makea  problem  (8)  nonlinear  and  one  may  apply  the 
concentratlon-coapactneaa  arguoMnta  to  thla  full  nonlinear  A-body  problem.  We  will  not 
puraue  thla  queatlon  here. 

We  now  conclude  this  section  by  examining  the  validity  of  HP  approximations  to  (8). 
In  (31]  I  E.  H.  Id.eb  and  B.  Simon  proved  (at  least  for  Coulos(blc  systems)  that  the  ground 
state  energy  given  by  (8)  and  Its  HP  approximation  E^p  given  by  (11)  have  similar 

asymptotic  behaviours  as  A  go  to  -h>.  Hore  precisely  one  haa  for  general  classes  of  V 


(62) 


E" 

.7/3 


HP 


E 


,7/3  A+»  TP 


A  A' 

where  is  the  Inflmum  of  the  so-called  Thomae-Permi  approximation  of  (8) 


(63) 


E  -  lnf{|i  /  p^^^dx  +  j  f  f  P(x)V(x-y)p(y)dxdy/ 

R  «R 

P  €  l’(*?)  n  L®^^(m^),  P  >  0  a.e.,  /  ^  dx  -  1} 


2  2/3 

where  y  is  given  by  (6i»  )  These  results  were  first  proved  by  E.  H.  Lleb  and 

B.  Simon  [31]  and  the  original  proof  was  later  simplified  by  B.  Baumgartner  [4], 

E.  H.  Lleb  (261,  (27]:  an  Inspection  of  the  proof  (confirmed  to  the  second  author  by 

E.  H.  Lleb)  shows  that  the  result  holds  for  general  V.  In  fact.  In  Huclear  Physics  It 
E 

HF 

is  expected  that  ~  <  0  (volume  energy  constant):  this  means  that  for 

realistic  V  E,pp  •  0.  Furthermore,  defining  E^p  by  (63)  where  /  .  pdx  •  1  is  replaced 

A  A  " 

,  E*^  ®HP'  ■'tF 

by  /  ,  pdx  “  A,  one  would  like  to  prove  that  — ’  —  —  C..  Of  course  If 

_3  A  A  A  0 

R 

E^  <  0,  then  one  deduces  from  (62)  that  E^/E^p  1. 

Concerning  the  TF  minimization  problems  (63),  let  us  mention  the  references  (36],  [33] 


where  related  problems  are  treated.  Applying  the  method  in  [33],  we  find  that  if  V  Is 


given  by  (5)  1^^  -  0  and  there  ie  no  minimum  while  if  V  is  given  by  (3)  or  (4)  the 

concentration-compactness  argument  applies.  And  it  is  shown  in  [33]  that  every  minimizing 
sequence  is  relatively  compact  up  to  a  translation  if  and  only  if  <0.  In  particular 

if  Eiff  <  0,  there  exists  a  minimum  and  if  a,  0  <  0  then  this  minimum  is  spherically 
syismetric,  nonincreastng  with  respect  to  |x|.  Finally,  one  checks  easily  that  if  a"  +  B' 
is  small  then  •  0  while  if  a,  B  are  negative  and  large  <  0. 


ri 


VIII.  Solutions  with  syametrias  of  Hartraa-Fock  aquations 


All  tha  minimisation  problams  wa  considarad  in  tha  pracadlng  sactions  ara  invariant 
undar  othoqonal  transformations  of  if  R  is  an  orthogonal  matrix  than  danoting  by 

V^(*)  *  S^CR*)  for  all  i  wa  chack  immadiataly  that 

for  all  tha  functionals  wa  considarad  pravioualy>  tdiila  tha  orthogonality  conditions  (10) 
still  hold  for  ). 

It  is  thus  natural  to  look  for  solutions  of  tha  HF  aquations  with  cartain  Invarianca 
properties  by  a  subgroup  of  the  group  of  orthogonal  transforms  of  R^.  For  Instance  one 
may  look  for  solutions  with  spherical  aymmatry  or  cylindrical  symmetry  and  in  particular 
one  may  study  the  same  HF  minimisation  problems  with  the  additional  constraint  of 
invariance  by  a  chosen  subgroup.  But  of  course  there  are  various  ways  to  impose  spherical 
symmetry  on  (or  cylindrical  symmetry).  One  possibility  is  to  Impose  that 

all  are  spherically  symmetrle  l.e.  really  depend  only  on  |x|.  However  this  is 

not  really  satisfactory  from  the  Physics  view  point  since,  even  when  solutions  with  such 
symmetries  exist  (and  this  is  not  always  the  case  in  view  of  numerical  experiments),  in 


general  such  a  solution  gives  a  value  to  tha  functional  which  is  too  high  to  yield  any 
information  on  problems  like  ( 8 ) . 

To  explain  the  meaning  of  spherically  symmetric  solutions  in  Nuclear  Physics  we  take 
an  example  namely  the  case  of  the  functional  (28)  and  to  simplify  we  assume  that  e  •  0  so 
we  consider 


=  . ’a’  “  T  -  2-  ^  ^  <>1^^  * 

-  J  ‘Vn  ^  Vp’  ^  I’Ppl'’  • 


Recall  that  we  work  with  the  following  othogonallty  conditions 

(65)  f  ,  »,s*dx  “5,,  for  1  <  i,  1  <  N  and  for  N+1  <  i,  j  <  A  . 


Now  If  w«  assume  there  exists  a  critical  point  of  B  given  by  (64)  with  the  constraints 

(65)  such  that  p^,  pp,  Tp  are  spherically  syometric  then  up  to  some  unitary  transform 
of  the  form  (29)  the  HF  equations  My  )»e  written  as 

(66)  -  dlv  ( — 5 - +  V  (r)e^  -  e^^Sj^  on 

2m  (r)  ’ 

<I 

where  e^,...,ej^  are  the  Imgrange  multipliers,  q-n  If  1<1<N,  q^p  If 

N-f1  <  1  <  X,  r  -  |x|  and  m^,  m^,  Vp  are  spherically  symmetric  functions  which  are 

easily  computed  from  the  expreeslon  of  B  (64).  The  quantities  ra^  are  often  called 

effective  masses.  It  Is  well-known  that  being  an  eigenfunction  of  the  elliptic 

^2 

operator  given  by  [“  dlv( — y—  V)  +  V  ]  lust  be  a  product  of  a  function  0Ar)  by 

2m^(r)  ’ 

a  spherical  harmonic  l.e.  an  eigenfunction  of  the  Laplace-Beltraml  operator  (-^g)  on 

the  sphere  S  of  R*'.  But  ttien  the  orthogonality  conditions  (65)  Imply  that  If  for  some 
1  €  (1,...,N}  (for  example)  and  “Ag'I'i  •  ®i’l'i  then  denoting  by  m  the 

multiplicity  of  the  eigenvalue  B^^  there  exist  (st-l)  Indices  In  distinct 

from  1  for  which  the  associated  ^  Is  also  an  eigenfunction  of  with  the 

eigenvalue  B^.  In  other  words  N  (and  Z)  splits  Into  the  sum  of  say  1(  multiplicities 
of  eigenvalues  of  (-dg)  and  the  angular  functions  4'^  associated  to  span  the 

elgenspaces  of  these  eigenvalues.  In  view  of  the  Increasing  multiplicity  of  eigenvalues 
of  (-&g)  as  they  increase,  It  Is  easy  to  see  that  for  given  N  and  z  they  are  only  a 
finite  nuBd>er  of  choices  for  the  angular  dependences  of  the  functions  This 

decong>o8ltlon  Is  precisely  the  meaning  of  a  spherically  symmetric  solution  for  HF 
equations. 

Of  course,  we  could  ailnlad.se  B  loosing  (65)  and  the  above  formulation  of  spherical 
symmetry  but  It  Is  somewhat  simpler  (and  better  for  the  values  If  E  we  get  this  way  a 
priori)  to  consider  instead 

I®  -  Inf{B[e^,...,»j^]  /  f «  h’(R^)  for  1  <  1  <  A,  (65)  holds, 

(67) 

/  jPTdx  <  •,  P|,'Pp''*n'^p  spherically  symmetric) 


Obaerv*  navarthaless  t..at  if  w«  find  a  minimum  of  (67)  than  by  tha  abova  arguaianta  tha 
alnlmum  la  raally  of  tha  form  wa  daacrlbad  abova  ao  thara  la  no  loaa  by  conaldarln^  tha 
mlnlmlalng  aat  daacrlbad  In  (67). 

hnd  one  provaa  eaally  the  following  raault  ualng  either  tha  concentratlon-coa^etnaaa 
argumanta  with  aymBMtry  ((  ] )  or  the  almplar  fact  that  If  p*  la  bounded  In  L^(>?)  n  H^(i^ 

and  la  apharlcally  aymaMtrlc  then  p*  la  coi^act  In  1^(1^ )  for  1  <  p  <  6  (aaa 
W.  Strauaa  [461.  H.  Baraatyckl  and  p.  L.  Llona  (5),  P.  L.  Llona  [37]).  Before  atatlng  the 
raault  wa  juat  naad  a  notation 

i®(m^,...,«j^)  -  inf{E[a^,...,aj^l  /  «  H^(«^),  /  ^ptdx  <  •  , 

R 

(68)  /  j  eor  1  <  1,  J  >  H,  H+1  <  1,  j  <  h  , 

P  ,  P  ,  T  ,  T  are  spherically  aymmatrlc) 
n  p  n  p 

where  >  0  for  all  1  e  ( 1 . . . . , A} . 

Theorem  VIII. It  Asaume  that  a  >  (6+i)/2,  a  +  SAy  >  S+0.  Then,  every  nlnlaialng  sequence 
of  (67)  la  relatively  compact  In  R^(R^)  (and  pr  Is  relatively  compact  In  1,^(R^))  If 
and  only  If  the  following  condition  holds 

(69)  I®  <  l®(m  ,...,m  ),  for  all  m  e  [0,1](1<1<A)  such  that  J  m,  <  A  . 
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In  particular  there  exists  a  minimum  of  (67)  If  (69)  holds. 

Remarks ;  t)  As  in  many  results  above,  the  condition  (69)  seems  difficult  to  check  for 
A  <  2  and  In  fact  numerical  computations  that  the  existence  of  a  minimum  Is  highly 
dependent  on  A. 

2)  Again  It  Is  possible  to  treat  the  case  when  h  -  0.  □ 

For  realistic  Interactions  and  HF  problems  (thus  Including  spln-orblt  forces)  the 
spherical  symmetry  Is  Imposed  by  considering  (v^,...,P^)  such  that  p^^,  Pp,  Tp  are 
spherically  symmetric;  Pj,(x,y),  Pp(x,y)  satisfy  Pq(Rx,Ry)  -  p(x,y)  for  all  rotations 
R  of  R^  and  for  all  q  -  n,p  and  J^,  Jp  have  the  form 

-r  |J„l<r)  , 


i 


M 


* 


and  with  thasa  oonstralnts  alnllar  raaulta  hold. 

Aa  wa  alraady  axplalnad  abova  two  arguatanta  «ay  ba  Invokad  to  prowa  tha  analoguaa  of 
Thaoraa  VXII.1  In  tha  oaaa  of  aora  raallatlc  Intaraetiona  T(x^).  ttthar  ona  appllaa  tha 
ganaral  arguaanta  of  P.  L.  Liona  [34]  ( eoncantratlon-ooapaotnass  prlnolpla  in  praaanea  of 
ayonatrlaa).  or  ona  aay  uaa  aora  atandard  coapaotnaas  arguaanta  d«a  to  apharleal  syaaatry 
aa  aantionad  abova.  In  tha  lattar  caaa  howavar  ona  naada  to  axplain  how  to  paaa  to  tha 
Halt  on  tha  tara 

/  /  V(x-y)|p  (x,y)|*d««r 

whara  la  tha  danalty  corraapondlng  to  a  alnlalalng  aaqoanoa  (f*,  ...,v*)  (thla 

boundad  In  No  thua  aaauaa  wa  hava  apharloally  ayaaMtrlc  oonflguratlona  l.a. 

p^(Rx,l^)  -  P,(x,y)  »  x,y  e  1^ 

for  all  notatlona  R  of  R^.  Slnca  V  dacaya  at  Infinity  It  la  anoagh  to  axplain  tthy, 

far  all  ((<*■,  |Pa(x>y)|^  ^|x-y|<R  coapaot  on  L^(R^xp^}.  wa  than  Introduoa 

C  (x)  -  /  lp„(x,y)|*dy 

“  |x-y|<R  * 

which  la  apharlcally  ayaaatrlc  (Cg,<Rx)  ~  Cg<x)).  And  la  boundad  In  W^'Vr^)  by  tha 
daflnltlon  of  P  and  bounda  on  a?.  By  p.  L.  Hona  [37],  wa  aaa  that  C_(x)  < 

X  *"  t  1 2 


on 


iPa  Ther«for«/  if  «•  prov«  that  la  boanAad  in  it  la  than  aaay  to 


coneluda  that  Cg,  la  compact  In  l4^(R^)< 

/  3  <  C  /  /  |v“(x)|^  |f“(y>|^dy/2ax 

R  R  1  lx-y|<R 

<  c  /  K  /  |Ri<x)|*  lf?(y)l*dy/*dx 

R  1  |x-y|<R 

<  C  I  /  3|a"(x)|C  /  |f“(y)I^dy/2ebc 

1  R  lx-y|<R 

<  C  I  (/  3l»®(x)l^dx)''2  (/  (  /  l»*(y)l*dy)dx)^2 

1  R^  ^  r'  |x-y|<R 

3 

<  C  A(4irR  )  2  , 


But 


for  varloua  constanta  C  >  0  and  wa  coneluda. 
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•  •j*-  -*• . 

-  W.  < 


"•"t  *  *V  *a*^  ■'.<  ■*> 


•  v. 


.  -wN' 


W«  now  conclud*  this  ssctlon  by  s  brlsf  dlsousslon  of  cyllndrlosl  syasstrio 
solutions.  Lst  us  dsnots  by  s  ■  >  s  «  X3  If  x  <■  (x^,X2/X3)  Is  s  gsnsxlo  point 

of  1^.  Argumsnts  similar  to  those  given  above  lead  to  the  following  problem 

1*  -  Inf{*[v^,...,»j^]/V^  c  v}{W?),  I  3  ptdx  <  (65)  holds, 

(70)  * 

P  ,P  are  functions  of  s,  s  only)  , 

n  p  n  p 

that  we  extend  In  the  following  class  of  problems 

I®(m^,...,m^)  -  Inf{K(e^,...,V|^l/f^  H^(«^)j  f  3  “  “l^lj 

(71)  for  1  <  1,  j  <  M,  M+1  <  1,  j  <  A»  /  3  ptdx  <  •»  Pn''’p'^n'% 

are  functions  of  s,  s  only)  , 

where  m^,...,mj^  are  nonnegative  constants. 

Again,  applying  the  coneentratlon-coi^ctness  arguments  we  see  that  If  a  >  {t*6)/i, 
a  -f  SAy  >  6'^6,  then  every  minimizing  sequence  of  (70)  Is  relatively  compact  up  to  a 
translation  In  z  If  and  only  If  the  following  condition  tMlds 

I*  <  i*(m^,...,mj^)  +  l*(1-m^,...,1-«j^)  for  all  m^  e  (0,1) 

(72)  ^ 

(1<1<A)  such  that  \  m.  e  (0,A)  . 
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and  there  exists  a  minimum  of  (70)  if  (72)  holds 


IX.  The  ihap*  of  th»  nttcl«u»  »nd  By— try  braltlnm. 

Admitting  that  tha  RF  approximation  la  valid«  than  tha  ground  atata  of  a  nuclaua  la 
auppoaad  to  ba  daacribad  by  tha  minima  of  tha  varioua  HF  minimization  probXama  studlad  in 
Sectlona  III  -  VI.  In  particular  tha  ahapa  of  tha  nuclaua  will  ba  datarmlned  by  tha 
danaity  pi  tha  nuclaua  ia  ^harical  if  P  ia  apharlcally  aymmatrlc,  or  mora  ganarally 
haa  tha  aynnatrlaa  that  p  poaaaaaaa. 

If  wa  heap  tha  notatlona  of  tha  pracading  aactlon,  wa  aea  that  the  spherical  symmetry 
is  broken  if  I*  >  1  but  it  nay  happen  that  1*  >  1°  •  I  in  which  case  the  spherical 
symmetry  is  broken  but  the  minimum  (if  it  exists)  still  preserves  the  cylindrical 
symmetry.  Hhlla  if  I*  >  I^  >  I  than  even  tha  cylindrical  symmetry  is  broken.  All  these 
phenomena  (and  many  others  related  to  s»re  elaborate  symmetries)  are  known  to  occur  in 
Nuclear  Physics  and  are  very  important,  flie  mechanism  Isehlnd  these  symmetry  breakings  is 
not  all  all  understood  neither  from  tha  Physics  viewpoint  nor  from  the  mathematical 
viewpoint.  He  propose  here  soma  vagus  explanations  on  soma  of  these  symmetry  breakings  and 
we  consider  as  examples  varioua  sluplar  modal  problaaw  which  could  help  understanding  these 
phenoBMna.  Before  going  into  these  exaivlas,  wa  would  Ilka  to  comment  on  the  physical 
meaning  and  isgilicationa  of  suoh  deformed  HF  ground  states.  Since  the  original  Hamiltonian 
is  rotationally  invariant  wa  know  that  tha  real  ground  state  should  have  a  "good  total 
angular  momentum".  This  aspect  whlcht  at  first  eight#  seems  to  be  a  defect  of  HF  theory 
hides#  on  the  contrary#  very  nice  features  aa  briefly  explained  in  the  following.  In  order 
to  restore  the  symmetry  the  HF  solution  is  raintarpretad  as  an  intrinsic  state  capable  of 
rotating  onto  intself.  Ihe  quantitation  of  such  collective  rotational  motion#  achieved  by 
A.  Bohr  [11]  a  long  tiSM  ago#  leads  to  a  oiodel  predicting  excited  states  whose  spectrum 

m2 

should  obey  the  simple  lawi  KI-M)  where  I  is  the  total  angular  momentum  and 

g  ia  the  Inartla  moswntum.  Such  typical  rotational  spectra  are  exhibited  by  serval  nuclei 
j152#154gj|^  154gj^  ^**Er)  and  tha  HF  theory  do  predict  in  their  cases  a  deformed  intrinsic 
structure.  It  must  be  pointed  out  hare  that  such  intarpretatlon  of  broken  symmetries  In 
terms  of  collective  swdes  is  currently  used  in  various  branches  of  modern  physics  (see  J. 
(lOldstone  modes  [22]#  Higgs  modes  in  non  Abelian  Gauge  theory  [24]).  Thus,  It  appears  that 


the  HF  method  le  a  much  more  powerful  tool  than  it  looks  a  priori  from  a  strict 
mathematical  vie%ipoint> 

We  begin  with  a  very  simple  example 

Example  1 ;  let  B  be  a  ball  centered  at  0  in  R^>  He  consider  the  minimisation  problem 
A  A 

I  ^  /g  «  I 


(73) 


for  all  1  <  1  <  A,  /g  u^Ujdx  ■  6^^  for  1  <  if  j  (  a} 


where  X  >  0  is  a  parameterf  u^(1<i<A)  are  real-valued  functionSf  f  is  a  continuous 
function  on  R**  satisfying  for  example 
(74) 


Urn  f(t)t”®'^®  <  0 


We  claim  that  if  A 


t++« 

2  or  if  A  »  3  then  for  X 


_  ill  the  density  p  «  J  uj  is  not 
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spherically  symmetric  where  (u^,...fUj^)  la  any  minimum  of  the  above  minimisation 
problem.  (The  existence  of  minima  is  a  standard  exercise  on  functional  analysis  since  we 
are  dealing  with  a  bounded  domain  B  and  the  nonlinearity  satisfies  some  appropriate 
growth  condition).  To  prove  this  claim,  we  denote  by  B(X)  the  value  of  the  above  infimum 
and  we  observe  that  for  X  -  0,  E(0)  is  nothing  but  the  sum  of  the  first  A  eigenvalues 
of  the  operator  -A  in  Hg(B)  and  that  the  corresponding  minima  for  A  -  2  or  3  are 
such  that  p  is  not  spherically  symmetric.  To  conclude'- we  just  have  to  prove  that  B(X) 
converges  to  E(0)  and  that  minima  (u^f...fuj^)  of  B(X)  converge  (extracting  enough 
subsequences)  to  the  minima  of  E(0).  Indeed,  observe  that 

E(X)  <  E(0)  4-01,  for  some  C  >  0 
while  (74)  implies  easily  that  for  some  C  >  0 

I  ^  fa  <  C 

i-1  ®  ^  ® 

X 


X  /„  f^p*)dx  0 


This  yields  on  one  hand  that  E(X)  ^4^  E(0)  and  on  the  other  hand  if  u^(1<i<A)  converge 
weakly  in  Hn(B)  to  u,  then 
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-r-  *rv** 


I'-Ji 


I  <  ii£  I  /bI^^iI^^ 

1-1  *  ^  vJoi-i  ®  ^ 


<  ii»  {  I  ♦  ^  /-  «"(P^)dx} 

X+0  1-1  ^  ■ 

<  11m  E(A)  -  «(0) 

X+0 

and  sine*  tiM  oonatralnta  pass  to  tha  lladt,  our  claim  la  provad. 


Ona  aaas  what  la  tha  machanlsm  Involvad  In  tha  abawa  ascam^la  and  It  saaaui  that  this 
machanlsm  plays  a  rola  In  Ihiclaar  Ihyslcst  roughly  spaaklng  It  la  axpactad  that  synawtry 
braaklngs  "hava  mora  chancas  to  occur"  for  thosa  h  such  that  tha  combinatorics  of  filling 
our  Slatar  datarmlnanta  with  apharlcal  harmonics  (as  wa  axplalnad  In  tha  pracadlng  sactlon) 
do  not  maka  poaslbla  tha  usa  of  only  tha  lo««ast  posslbla  alganvaluas  (or  anargy  lavals) 
of  (-A^).  This  axplanation  la  vary  snich  ralatad  to  what  Is  callad  In  Rtyalcs  suiglcal 
nuB^rs«  Of  couTSSf  this  tantatlva  axplanation  has  to  ba  conflrmad  or  Inflrmad  by  tha 
axamlnation  of  sera  raallstlc  problams  than  (73).  Wa  proposa  anothar  nodal  problam  for 
which  It  would  alraady  ba  Intarastlng  to  daclda  Whathar  thara  Is  symmatry  braaklng  or 
not.  Wo  will  only  mantlon  tha  casa  A  -  3. 

Bxawpla  2i  Wo  consldar  now 


1(1,1)  «  Inf{/  ,|Vu|^  +  jVvj^dx  -  /  ,f  u*(x)V( x-y)v*(y)dxdy  / 


u,  V  «!(’(■*),/  ,  u  dx  -  /  ,  v*dx  -  1} 

vr 


whara  X  >  0  and  V  Is  spharically  symmatrlc  and  satlsflas  V  -  V2  whara 

V  e  L^^(a^)(l  -  1,2)  tor  soma  €  [-j,  -) .  iha  main  dlffaronca  batwaan  (75)  and  HF 
problams  Is  tha  fact  that  wa  do  not  assumo  anynora  that  j ^  uv  dx  -  0  and  It  Is  posslbla 
tliat  this  typo  of  constraints  plays  an  Important  rola  In  symmatry  braaklngs.  Wa  provo 


balow  that  as  soon  as  1(1,1)  <  0  all  minimizing  saquoncas  ara  ralatlvoly  compact  in 
H^(B^)  up  to  a  translation  and  thus  thara  exists  a  ad.nlmum  of  (75).  By  symmatrlzatlon 


•:\  'v;v;v  v  ^ 


.a.  a.  ^ 


t«chiiiqiMS  mlnlaa  ar«  apherically  synmtric  if  V  ia  nonnagatlvaii  nonlneraaaing  with 
raapact  to  |x|.  The  caaa  whan  V  doaa  not  hava  thaaa  propartiaa  la  totally  opan. 

To  prova  tha  abova  clalur  wa  hava  to  ahow  (ualng  tha  coneantratlon-compactnaaa  nathod) 
that  If  1(1,1)  <  0  than 

1(1,1)  <  1(1, u)  >  1(1-1, 1-u)  for  0  <  1  <  1,  0  <  M  <  1  and  0  <  X*u  <  2 
whara  1(1, y)  atanda  for  tha  aaaM  Inflnua  aa  In  (75)  but  with  tha  conatralnta 

/  u^dx  -  1,  /  3  v^dx  -  u.  Iha  proof  of  t)iaaa  atrlct  Inaqualltlaa  uaaa  tha  fact  that  If 

«  * 

1-0  or  If  u  -  0  than  1(1, w)  -  0  whlla  If  0<1<1,  0<tJ<1  than  1(1, u)  ~  Iv 
B({,  -i)  where 

K(a,t)  -  Inf{a  /  |7u|^dx  +  t  /  [Vul^dx  -  /  /  u^(x)V(x-y)v*(y)dxdy  / 

U,  V  «  H^m^),  /  tt^dx  -  /  ,  v*dx  -  1} 

«  » 

and  thua  E(8,t)  la  nondecreaalng  with  reapact  to  a  or  to  t>  Cbaervlng  naxt  that  ly  * 
(1-1) (1-U)  <1  If  0<l<1,0<y<1  wa  deduce  the  above  atrlct  aubaddltlvlty 
Inequality  by  renarlclng  that 

8(1,1)  <  l(p  -J-),  *(1,1)  <  *(—,  ~)  . 

□ 

We  conclude  thla  aectlon  by  nentlonlng  that  the  atudy  of  varloua  nuclell  aeaaia  to 

Indicate  that  the  nechanlau  wa  llluatrated  by  the  alaple  example  1  apparently  doaa  not 

cover  all  the  poaalble  waya  the  apherlcal  ayoBMtry  la  brolcen.  To  explain  thla  claim,  let 

ua  flrat  explain  how  apln  dependence  (and  apln-orblt  forcea)  makea  the  above  deacrlptlon  a 

bit  more  complicated.  Indeed,  In  auch  a  caae,  the  aequence  of  HF  levels  la  typically  of 

the  following  form  (I8V2  multiplicity  2,  Ip  3/2  multiplicity  4,  Ip  V2  multiplicity  2, 

Id  5/2  multiplicity  6,  2a  V2  multiplicity  2...)  where  the  atatea  (levela)  are  labelled  aa 

it  la  uaual  by  the  set  n,  t,  j  where  j  denotes  the  eigenvalue  of  the  total  angular 
♦  +  ♦ 

momentum  j  «  e  s.  Thua  the  degeneracy  (multiplicity)  of  the  level  (n,i,j)  Is  2ji-1. 

12  2B 

With  this  scheme  one  checks  that  both  nuclei  ,  C  and  .  .8  have  nucleon  numbers 

o  14 


coaq^tlble  with  a  spharlcal  HF  solution  built  on  ths  lowsst  levels  or  eigenvalues.  Yet, 

the  HF  cosq^tations  lead  to  solutions  which  are  not  spherical  but  defomed.  In  fact  they 

correspond  to  axially  defomed  shapes  i.e.  solutions  with  cylindrical  syoanetry 

2^2  ,2 

(X  4y  2  ^  c  c 

— r*—  “J  ”  1 J  which  are  prolate  (I.e.  T  >  oblate  corresponding  to  —  <  1)  -  see  M 

a^  c^  * 

Glrod  and  B.  GreuMiatlcos  [21]}  K.  Kuaar,  Ch.  Lagrange,  M.  Glrod  and  B.  Gramoatlcos  [25]. 


X.  External  field  method 


To  present  the  external  field  method  it  is  trarth  saying  a  few  words  on  the  numerical 
computation  of  HF  minimization  problems.  Because  of  the  quite  complicated  form  of  the 
functionals  one  has  to  minimize,  some  numerical  methods  which  are  currently  used  (typically 
Galerkin  type  methods  based  on  spherical  harmonics,  or  two  basis  of  spherical  haxvonics 
centered  at  different  points...)  make  difficult  to  avoid  symmetries  and  seem  to  favor  the 
possible  local  minima  with  spherical  or  cylindrical  symmetry.  And  in  all  cases  the 
numerical  methods  break  the  translation  invariance.  These  remarks  explain  why  the  standard 
problem  in  the  minimization  of  nonconvex  functions  of  avoiding  local  minima  in  order  to 
find  the  absolute  minimum  seems  even  more  acute  in  KF  problems.  One  idea  to  avoid  this 
difficulty  is  to  deform  the  shape  of  the  density  by  an  external  field  acting  on  the  system 
as  an  additional  constraint.  As  we  will  explain  below  this  approach  is  not  only  useful  for 
numerical  purposes  but  is  also  relevant  for  physics. 

To  explain  the  principle  of  the  external  field  method,  we  consider  a  functional 

on  a  manifold  N,  bounded  from  below  and  we  are  interested  in  the  minimization  problem 
(76)  E  -  Inf(E(u)|u  €  M}  . 

Now  if  Q  is  some  given  functional,  and  q  e  R  we  consider  the  same  minimization 

problem  where  we  add  the  constraint  Q(u)  ==  q  (which  in  HF  problems  represents  the  action 
of  the  external  field) 


(77) 


E(q)  =  Inf{E(u)  /  u  €  M,  Q(u)  “  q}  . 


It  is  obvious  that  E  >  Inf  E(q)  and  that  if  this  infimum  is  achieved  at  qQ  and  if 

qeR 

E(qg)  is  achieved  at  some  Ug  then  Ug  is  a  minimum  of  (76).  But  observe  also  that  if 


we  assume  that  at  some  qg  E  is  differentiable,  E'(qg)  «  0  and  that  for  q  near  qg 


E(q)  is  achieved  at  some  u^^  differentiable  with  respect  to  q  at  qg,  then  u^  1 


critical  point  of  E.  Indeed,  u  being  a  minimum  of  E(q. ),  there  exists  a  Lagrange 

40  '' 


multiplier  9  such  that 
(78) 


E-(u,^)  -  9Q'(u,^) 


du 


while  if  we  denote  by  v_  = 


then  differentiating  the  relations  Q(u  )  **  q. 
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XI.  Ti— -d«p«nd«nt  Hartree-Fock  xiuationa. 

ni«  tim*-d«p6nd«nt  Hartrea-Fock  aquationa  (TDHF  tn  short)  ara  couplad  nonllnaar 
SchrSdlngar  aquations.  Glvan  any  HF  functional  B(a^/. . . /Sj^l  as  In  tha  pracadlng 
sactlons>  TDHF  aquations  nay  ba  wrlttan  as  follows 


(82)  ^  3t^  ■*■  If"  -  0  on  W?  »  (0,-),  1  <  k  <  A  . 

Of  coursa,  to  solva  (82)  ona  has  to  add  Initial  conditions 

(83)  ay(x,0)  -  a^(x)  on  I^,  1  <  k  <  X 

whara  ara  glvan  ( 1<k<A) . 


For  axampla  If  X  la  glvan  by  (21),  than  (82)  may  ba  rawrlttan  as 

®*k  fl  3 

i  3^ - m"  ^’k  “  7  +  "•x  -  on  R  >«  (0,-),  1  <  k  <  A 

d  Y  2 

whara  W  •  -up  +  'jT-'jip  +  6p  . 

For  tha  motivations  In  Huclaar  Physics  for  studying  TDHF,  wa  rafar  tha  raadar  to 
H.  Flocard  [16]. 

Wa  will  not  give  rasults  concarnlng  the  resolution  of  the  Cauchy  problem  (82),  (83)t 

let  us  just  mention  tha  trarks  by  J.  Glnlbra  and  G.  Valo  [18],  [19],  [20] >  for  all 

Interactions  except  Skyrrae's  there  la  no  special  difficulty  to  aolva  (82)  -  (83).  Many 

mathasMtlcal  rasults  on  systems  Ilka  (82)  ara  based  on  tha  various  conservation  laws 

satisfied  by  solutions  of  (82):  for  example  multiplying  (82)  respectively  by  and 

®’k 

taking  the  Imaginary  part  and  by  and  taking  tha  real  part  ona  finds  Integrating 

over  R^ 

(84)  ^  '  1<k<A 

(85)  E[a^,...,a^]  Is  Independent  of  t  . 

,  * 

Similarly,  one  obtains  that  J  ,  a  g.dx  Is  Independent  of  t  for  1  <  k,  t  <  A. 

r3  K  l 

He  next  observe  that  solutions  of  HF  aquations  (up  to  unitary  transform)  lead  to 
stationary  solutions  of  TDHF  aquations  where  stationary  means  that  p,  t  are  Independent 


of  t:  more  precisely  we  have  seen  In  the  preceding  sections  that  we  may  write  the  HF 


equations  as 


. ’a’  ■  Vk  ‘  "  '■  • 

Ityt 

Then  obviously  -  e  <  1<k<A)  defines  a  solution  of  (82). 

In  particular,  any  minimum  of  the  HF  minimization  problems  leads  to  a  stationary 
solution  of  TOHF  equations.  It  is  sho«m  in  T.  Cazenave  and  P.  I.  Lions  [12}  that,  if  the 
subadditivity  conditions  given  in  the  preceding  sections  via  the  concentration-compactness 
arguments  hold,  minima  of  HF  minimization  problems  are  orbltally  stable  in  TDHF 
equations.  Let  us  also  point  out  that  similar  arguments  show  that  minima  of  the  HF 
minimization  problems  with  additional  symmetry  constraints  (see  section  VIII)  are  orbitaliy 
stable  with  respect  to  perturbations  with  the  same  symmetries. 

But  since  all  solutions  of  HF  equations  lead  to  stationary  solutions  of  TDHF 
equations,  the  study  of  all  possible  solutions  of  HF  equations  presents  some  interest.  In 
particular  one  may  look  for  critical  points  of  E(e^, . . .  ,9j^)  with  the  additional 
orthogonality  constraints.  The  only  approach  we  know  one  might  try  is  through  min-max 
principles  as  it  is  done  in  H.  Berestycki  and  P.  L.  Lions  [51,  P.  L.  Lions  [38],  [39]  for 
related  problems.  This  approach  requires  spherical  symmetry  of  the  functions  (p^,...,e^). 
We  only  have  very  partial  existence  results  in  that  direction. 

In  fact  Nuclear  Physics  considerations  (collisions  of  heavy  ions)  indicate  that  it 
would  be  interesting  to  find  all  periodic  solutions  of  TDHF;  again  one  has  to  define  the 
precise  meaning  of  periodic  solutions.  For  example,  if  no  differences  are  made  between 
neutrons  and  protons,  a  solution  of  (82)  is  said  to  be  periodic  of  period  T  if  there 
exists  a  unitary  transform  U  such  that 

(»,(T),...,9j(T))  =  0(9,(0),...,<p^(0))  on  . 

Observe  that  this  implies  that  the  densities  p,  t  are  indeed  periodic  of  period  T  (in 
the  usual  sense).  In  fact,  an  even  more  general  (possibly)  notion  of  periodic  solutions 
consists  in  requiring  the  density  p(x,y)  to  be  periodic.  It  seems,  at  least  numerically, 
that  there  are  many  periodic  solutions  of  TDHF  eqpiatlons  and  this  is  another  major  open 


question 


A  final  remark  on  this  topic  concerns  the  possibility  of  having  stationary  solutions 
of  TDHF  equations  which  are  not  obtained  through  solutions  of  HF  equations.  We  Illustrate 
this  possibility  on  a  simple  example  of  a  system  of  two  nonlinear  equations. 

Example;  We  consloer  the  following  system  of  two  coupled  nonlinear  Schrbdlnger  equations 


lip^  -  Ae  =  p 


Y-1 


on  R  *  ( 0 ,» ) 


(86) 


li(i^  -  Ai(i  =  t|i  on  R^  X  (0,«>) 


where  p  «  (|»P  +  and  1  <  y  <  5/3. 

Let  (1),  m  >  0>  we  look  for  solutions  of  (86)  of  the  form 
^  9(x,t)  ”  i^“^(coa  mt  u(x)  +  sin  mt  v(x)) 

I  i|i(x,t)  =  e^“'^(-8ln  mt  u(x)  +  eos  mt  v(x)) 


And  we  find  the  following  nonlinear  system  for  u,  v 

Y— 1  3 

-Au  +  Imv  +  oiu  -  p  u  In  R 

(87) 


Y-1  3 

-Av  -  Imu  +  uv  «  p  V  In  R  . 

Observe  that  If  we  show  there  exist  solutions  (u.v)  of  (87)  then  the  aliove  9(x,t), 
i|)(x.t)  yield  stationary  solutions  of  (86)  (p  Is  Independent  of  t).  However,  If  9 ,  i|> 
were  solutions  of  (86)  built  through  the  stationary  analogue  of  HF  equations,  this  would 
Imply  that  v  «  -lu.  Thus,  we  want  to  exhibit  solutions  of  (87)  with  v  ?  -  lu.  To  this 
end  we  consider  the  following  minimization  problem 


(88) 


I  =  Hin{/  ,|Vu|^  +  |Vv|  +  w|u|  +  “|v|  +  Re(im  vu)dx/ 


u,  V  e  H  (R  ),  /  p^dx  =  1} 


This  problem  is  solved  as  in  [33]  by  the  concentration-compactness  arguments  and  thus 
there  exists  a  minimum  (u,v)  of  (88).  Now,  if  we  have  v  =  -lu,  this  would  Imply 

I  >  Min{/  2|Vu|^  +  2u|u|^  +  m]  u|  ^dx/ueH^  (R^  ) ,  /  |u|^dx  =  2*^} 

R^  R^ 

and  this  last  minimum  is  strictly  larger  than 
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In  this  aection,  we  want  to  praaent  a  dlffarant  approach  to  the  atudy  of  nuclell 
naaely  the  ao-callad  Hartree-Fock-Bogolyubov  approximation  (M.  N.  Bobolyubov  [10], 

C.  Bloch  and  A.  Naaalah  [9],  J.  B.  Bardeen,  L.  H.  Cooper  and  J.  R.  Schrleffer  [3]).  We 
will  not  attempt  here  to  explain  the  Hartree-Fock-Bogolyubov  method  and  we  refer  to  the 
intereated  reader  to  J.  Decharge  and  D.  Gogny  [14],  j.  G>  Valatln  [47].  We  will  only 
deacrlbe  the  typical  adnledaatlon  problem  arlalng  In  thla  theory  which  may  be  thought  of  aa 
an  Improved  approxlaiatlon  of  the  A-body  problem  conaldered  In  aectlon  II.  To  almpllfy,  we 
will  again  Ignore  the  apln  dependence  and  apln-orblt  forcea.  The  minimizing  aet  la  given 
by 


M  •  /  J  u^v  +  u.Vj^dx  -  0 

/  ,  u  u*  +  V  v*dx  «  6  for  all  l,j  >  1»  I  /  |v  |^dx  «  A>  . 
R  J  J  ’  1>1  R 


We  now  Introduce  the  Hartree-Fock-Bobolyubov  problem  (HFB  In  abort) 

^  I  7  I  f  P(x)V(x-y)p(y)dxdy  + 

^  1>1  ^  *  .3^.3 

R  x|t 

-  4-  /  /  V(x-y)|p(x,y)|^dx  +  i  /  /  V(x-y)  |  K(x,y)  |  ^dxdy 

_3^_3 

R  *R  R  *R 

where  P(x)  -  I  |v  (x)|^,  P(x,y)  -  J  v  (x>v*(y),  ic(x,y)  -  J  u,(x)v*(y)  for  x,  y  e  iP 
1>1  ^  1>1  i>1 

And  one  may  chooae  for  example  (aa  In  [14])  the  potential 

V(x)  -  I  Wj  exp(-|x|*/uj) 

1-1,2 

where  are  conatanta. 

Moat  of  the  reaulta,  remarka  and  open  problema  given  In  the  preceding  aectlona  may  be 
adapted  to  the  study  of  problem  (89).  Before  explaining  how  we  may  apply  the 


concentration-compactness  program  on  this  problem,  let  us  first  mention  the  Important 


II: -.1  ".-I  «.TW 


K 


connection  between  the  above  problem  and  more  atandard  HF  problems*  Deriving  the  Euler- 
I*agrange  equations  of  the  above  minimisation  problem  (which  by  the  way  are  called  Hartree- 
Fock-Bogolyubov  equations)  one  sees  immediately  that  if  •  •  ,Vj^)  is  a  solution  of  the 

HF  equations  then  choosing  Vj^  =  v^^  for  1  <  1  <  A,  Vj^  5  0  for  i  >  A,  Uj^  s  o  for 
1  <  1  A,  Uj^  arbitrary  satisfying 

/^3  u^«*dx  - 

for  1  +  A  <  i,j,  we  find  a  particular  ("trivial")  solution  the  HFB 

equations.  In  soma  vague  sense,  the  HFB  problem  contains  the  HF  problem. 

We  now  conclude  with  a  brief  explanation  on  the  way  we  may  apply  the  concentration- 
compactness  arguments  to  the  above  problem.  We  apply  the  usual  concentration-compactness 
lemma  to  the  density  p.  And  we  see  that  minimising  sequences  the  above 

minimization  problem  are  compact  up  a  translation  in  Ii^(R^)  x  H^(R?)  if  and  only  if 

1  <  1(M^,M2,W)  +  I(M5,HJ,-N) 

for  all  hermltlan  matrices  H^,  Mj#  MJ,  N  satisfying  Mj,  >  0, 

Tr(M^)  «  [0,A],  Tr(M')  -  A  -  Tr(M,),  M,  +  Mj  +  M*  +  MJ  -  1  and  |  (IJ5 ,0  I  <  (M,C  ,5  )’'^2 
(H2C,5)*''^2  ,  w  C  €  where  in  addition  0  »<  M,  +  M2,  M,  +  1*2  ^  the  above  inequality 

X  •  Xnf {x (N^ ,M2 , 0 )  /  ^  Of  M^  •  M^,  M2  •  M2,  M^  +  M2  •  if  Tr(M|)  ■  A}  and  finally, 

denoting  by  F(l,j)  the  1-j  component  of  the  matrix  F  (1<l,j),  the  definition  of 
I(M^,M2,W)  Is  given  by 


I(M,,M2,N)  -  lnf{/  3^  iFv^l^dx  + -J  /  /  P(x)V(x-y)p(y)dxdy  + 

*  R^xr’ 

-w  S  !  V(x-y)  |p(x,y)  l^dxdy  + -1  /  /  V(x-y)  |  )i(x,y)  |  ^dxdy  /  (u”,v  )  t  x  h\ 

•.exa  *<»•»  ** 


-3,  =3 

R  XR 


R^xR^ 


/  3  Fj^v*dx  -  M^(i,J),  /  3  Uj^u^dx  -  M2(i,j),  f  3  Uj^v*  +  -  M(irj)}  • 

In  fact,  as  in  the  case  of  HF  problems,  using  the  various  invariances  of  the  above  problems 
it  is  possible  to  restrict  the  above  strict  subaddltivlty  inequalities  to  a  particular 
class  but  we  will  not  pursue  this  matter  here. 
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